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19.  ABSTRACT  (con't) 


examinations  were  made  of  both  the  as-recieved  DU  and  DU  specimens  subjected 
to  large  dynamic  shear  deformations. 

The  contained  fragmenting  cylinder  (CFC)  technique  for  determining  shear 
band  kinetics  parameters  is  reviewed  and  the  shear  band  size  distributions 
obtained  from  CFC  experiments  on  both  RHA  and  DU  are  presented.  A  preliminary 
set  of  RHA  shear  band  kinetics  parameters  and  a  revised  set  of  4340  steel 
parameters  were  obtained  for  the  modified  SHEAR3  computational  shear  band  model. 

The  development  and  testing  of  a  two-dimensional  wave  propagation  code 
(C-HEMP)  particularly  well-suited  to  armor  penetration  simulations  has  been 
largely  completed.  C-HEMP  currently  incorporates  a  complex  slide  line  capability 
(for  non- intersecting  slide  lines),  a  rezoner  that  handles  large  deformations, 
and  the  SHEAR3  model.  C-HEMP  with  SHEAR3  was  used  to  simulate  previously  performed 
armor  penetration  experiments,  and  C-HEMP  correctly  predicted  the  penetration 
and  fragmentation  phenomenology  observed  at  impact  velocities  near  and  far  above 
ballsistic  limit,  although  excessive  cell  distortion  currently  prevents  complete 
simulations  of  deep  penetration  cases. 

Details  of  the  SHEAR3  model  are  presented  in  Volume  2. 
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EXECUTIVE  SUMMARY 


This  is  the  first  volume  of  a  two-volume  series  comprising  the 

third  annual  report  for  the  program  "Computational  Model  for  Armor 

5>- 

Penetration,  performed  by  SRI  International  for  the  Ballistic  Research 
Laboratory  (Contract  No.  DAAK11-78-C-0115) .  ^The  objective  of  the  pro¬ 
gram  was  to  develop  a  phenomenologically  sound  material  disintegration 
model  for  computationally  simulating  armor  penetration  by  projectile 
impact.  The  model  would  provide  the  capability  to  predict  the  ballistic 
limit  and  downrange  fragment  size  and  velocity  distributions  and  would 
eventually  be  used  by  the  DoD  for  armor  and  penetrator  design.  The 
materials  studied  in  the  program  were  4340  steel^-fftc^)  anc*  depleted 
uranium  alloy,  (DU)  for  the  penetrators,  and  rolled  homogeneous  armor 

(RHA)  for  the  target. 

/  > 

—la  Volume  ~1  we  reports  on  the  experimental  and  analytical  work 

b 

performed  during  the  third  year  Co  conclude  the  development  of  a 
computational  model,  SHEAR3,  that  describes  the  dominant  material 
failure  and  fragmentation  process  (adiabatic  shear  banding)  responsible 
for  long-rod  penetrator  erosion  and  thick  armor  plate  plugging  and 
fragmentation.  Volume  1  also  describes  application  of  SHEAR3  with  the 
Lagrangian  wave  propagation  code  C-HEMP  to  computationally  simulate 
penetration  of  armor  plate  by  normally  impacting  long  rods  both  near  the 
ballistic  limit  and  at  higher  velocities  sufficient  to  produce 

f  I  ■'  j 

significant  back-of-the-armor  fragmentation,  i  h. 

Volume  2  documents  in  detail  the  SHEAR3  model. 


b 
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In  the  following  paragraphs  we  summarize  both  the  third  year 
results  and  the  overall  results  of  the  three-year  program. 

Our  effort  to  obtain  material  property  data  for  calibrating  the 
SHEAR3  model  was  completed  this  past  year  with  the  f ractographic 
examinations,  and  contained  fragmenting  cylinder  (CFC)  experiments  on 
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DU,  and  with  the  quantitative  analysis  of  the  CFC  experiments  on  both  DU 
and  RHA  to  obtain  shear  band  size  distributions.  Using  these  and 
previousy  determined  results  for  4340  steel,  we  performed  a  series  of 
computational  simulations  of  CFC  experiments  to  obtain  improved  esti¬ 
mates  of  the  SHEAR3  parameters  for  both  RHA  and  4340  steel. 

The  SHEAR3  subroutine  must  be  used  in  a  wave  propagation  hydrocode, 
and  we  have  developed  a  code,  C-HEMP,  that  is  particularly  well-suited 
for  armor  penetration  calculations  in  which  sophisticated  failure  models 
such  as  SHEAR3,  must  be  used.  During  the  third  year  of  the  program, 
C-!'.IMP  was  debugged,  tested,  and  partially  documented  (see  Volume  3). 

Finally,  SHEAR3  was  used  in  C-HEMP  to  computationally  simulate  the 
normal  impact  of  long  4340  steel  rods  into  RHA  plates  at  velocities 
ranging  from  near  to  well  above  the  ballistic  limit. 

The  status,  after  three  years,  of  our  attempt  to  understand  and 
predict,  in  detail,  armor  penetration  and  disintegration  based  only  on 
known  material  properties  is  as  follows: 

The  observed  behavior  shows  what  appears,  at  first  glance,  to  be  an 
enormous  variety  of  effects  ranging  from  simple  cratering  and  back  sur¬ 
face  bulging  to  plugging  (see  Figure  34,  Volume  1)  to  penetration  with 
catastrophic  back-of-the-armor  fragment  sprays  (see  Figure  33,  Volume 
1).  However,  the  first  significant  result  of  the  program,  reported  in 
the  first  annual  report,  was  that  all  the  observed  significant  material 
failure  and  disintegration  (erosion  of  the  penetrator  nose,  plugging  of 
the  armor,  and  downrange  fragmentation)  were  due  to  a  single  material 
failure  mode,  namely  adiabatic  shear  banding.  Thus,  our  modeling  effort 
centered  almost  immediately  on  developing  a  model,  SHEARS,  that  is  based 
on  measurable  laboratory  material  properties  and  that  describes  the 
nucleation,  growth,  and  coalescence  of  adiabatic  shear  bands  to  form 
either  one  large  fragment  (the  plug)  or  many  small  fragments  (penetrator 
nose  disintegration  and  erosion,  as  well  as  back-of-the-armor 
fragmentation) . 


For  normal  impact  SHEAR3  appears  to  correctly  predict  penetrator 
nose  erosion,  the  onset  of  plugging,  and  basic  features  of  high  velocity 
fragmentation  correctly.  However,  excessive  cell  distortion  prevents 
C-HEMP  from  carrying  the  calculation  to  completion  for  thick  armor 
plates.  Nonetheless,  the  calculated  size  distribution  of  the  fragments 
surrounding  the  nose  of  the  penetrator  when  the  calculation  halts  is 
very  similar  to  that  observed  experimentally  for  the  downrange  fragment¬ 
ation,  as  reported  in  our  second  annual  report. 

The  calculations  have  added  significantly  to  our  understanding  of 
the  penetration  process,  particularly  the  nature  of  the  transitions  from 
armor  bulging  to  plugging  to  disintegration.  These  transitions  proceed 
as  follows. 

At  velocities  approaching  the  ballistic  limit,  regions  of  complete¬ 
ly  shear  banded  material  coalesce  to  form  lobes  that  extend  towards  the 
rear  of  the  armor  plate;  eventually,  this  material  breaks  out  to  form  a 
plug  (See  Figures  38  and  37,  Volume  1).  A  key  feature  of  this  process 
is  that  the  impact  velocity  is  below  that  of  the  propagation  velocity 
for  significant  plastic  strain  levels.  Thus,  the  back  of  the  armor 
bulges  as  information  about  the  impact  reaches  it.  If  shear  banding 
were  absent,  this  back  surface  bulge  would  be  relatively  smooth,  but  the 
shear  banding  causes  localization  of  the  strain  at  the  incipient  plug 
boundaries  (See  Figure  40,  Volume  1).  Thus,  the  ballistic  limit  is 
determined  by  the  transition  from  bulging  to  plugging,  and  is  caused  by 
shear  banding  in  lobes  whose  position  is  affected  by  the  penetrator 's 
communication  with  the  rear  of  the  armor  plate.  Simultaneously,  the 
shear  banding  determines  the  penetrator  nose  erosion  (see  Figure  37, 
Volume  1). 

At  velocities  well  above  the  ballistic  limit,  another  transition, 
occurs:  from  plugging  to  significant  back-of-the-armor  fragmentation. 

When  the  impact  velocity  exceeds  the  propagation  velocity  of  significant 
plastic  strain,  the  penetrator  can  no  longer  communicate  with  the  rear 
of  the  armor  plate,  and  much  less  bulging  occurs.  Thus,  the  large 
plastic  deformation  and  associated  material  disintegration  from 
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coalesced  shear  bands  are  confined  within  a  "bow  wave"  around  the 
eroding  penetrator  nose  (see  Figure  36  and  compare  with  Figure  37). 

When  this  bow  wave  reaches  the  rear  of  the  plate,  a  spray  of  fragments 
is  released. 

In  summary,  SHEAR3  in  C-HEMP  has  significantly  clarified  our  under¬ 
standing  of  the  penetration  process  for  long  rods  impacting  thick  armor, 
and  has  given  us  the  capability  of  studying  the  sensitivity  of  the 
ballistic  limit  and  downrange  fragmentation  to  changes  in  material 
properties.  However,  full  use  of  this  capability  must  await  further 
development  of  C-HEMP  to  allow  computational  simulations  of  oblique 
impacts  and  of  deeper  penetrations.  Recommendations  for  such  further 
development  are  given  in  Chapter  VI  of  Volume  1. 
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PREFACE 


This  report  is  Volume  1  of  a  two-volume  series  representing  the 
third  annual  report  on  the  program  "Computational  Model  for  Armor 
Penetration,"  performed  by  SRI  International  for  Ballistic  Research 
Laboratory  (Contract  No.  DAAK11-78-C-0115) .  The  titles  and  authors  of 
the  individual  volumes  in  this  series,  as  well  as  those  of  the  previous 
annual  reports  for  this  program,  are  as  follows: 

_ Title _  _ Authors _ 

Vol.  1  -  Computational  Model  for  D.  C.  Erlich,  L.  Seaman, 

Armor  Penetration  T.  Cooper,  D.  R.  Curran, 

and  R.  D.  Caligiuri 

Vol.  2  -  Development  of  a  Model  for  L.  Seaman 

Shear  Banding:  SHEAR3 

First  Annual  Report  - 

Computational  Model  for  D.  R.  Curran,  R.  Burback, 

Armor  Penetration  R.  D.  Caligiuri, 

D.  C.  Erlich,  L.  Seaman, 
and  D.  A.  Shockey 

Second  Annual  Report  - 

Computational  Model  for  D.  C.  Erlich,  L.  Seaman, 

Armor  Penetration  R.  D.  Caligiuri,  and 

D.  R.  Curran 

This  program  has  benefited  greatly  from  the  technical  support  of 
Dr.  Gerald  Moss,  the  BRL  project  monitor,  and  Drs.  Tim  Wright  and 
Richard  Vitall  of  BRL. 

The  authors  wish  to  acknowledge  the  excellent  experimental  and 
computational  support  at  SRI  International  provided  by  A.  Urweider, 

D.  Henley,  D.  Petro,  H.  Hanna,  B.  Lew,  and  J.  Kerapf. 
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I  INTRODUCTION 


The  traditional  empirical  approach  for  designing  armor  and  pene- 
trator  systems  has  often  been  effective  in  gaining  successive  incre¬ 
mental  improvements  over  the  established  state  of  the  art.  However,  as 
more  complex  systems  are  introduced,  the  very  lengthy  and  expensive 
empirical  approach  becomes  less  and  less  tractable.  Recently  developed 
computational  methods  can  now  handle  complex  shapes,  deformations,  and 
material  behavior  and  are  therefore  more  likely  to  be  effective  in 
streamlining  armor  and  penetrator  design. 

The  computational  approach  to  armor  and  penetrator  design  is  an 
iterative  procedure  beginning  with  approximate  computer  simulations  of  a 
penetration  scenario,  using  existing  simple  models  for  material  behavior 
and  estimated  values  of  material  properties.  Test  firings  are  then  per¬ 
formed  to  evaluate  the  performance  of  the  design,  to  provide  data  for 
comparison  with  the  calculated  results,  and  to  indicate  the  adequacy  of 
the  material  behavior  models  and  material  property  values.  The  test 
results  suggest  changes  in  the  design  and  guide  the  selection  and  devel¬ 
opment  of  dynamic  material  property  tests  and  improved  material  behavior 
models  for  the  next  series  of  computer  simulations.  Test  firings  of  a 
second  generation  design  are  then  performed,  and  the  results  are  com¬ 
pared  with  the  computational  results,  as  before. 

The  procedure  is  repeated  until  the  system  attains  the  desired  per¬ 
formance.  This  approach,  currently  being  developed  by  DoD,  reduces  the 
overall  design  effort  because  the  interaction  of  experiments  and  calcu¬ 
lations  lead  more  rapidly  to  an  acceptable  design  than  a  purely 
empirical  approach. 

Computer  simulation  of  armor  penetration  must  model  the  following 
material  behavior:  a  thermodynamic  equation  of  state,  the  plastic  flow, 
and  the  material  disintegration.  Existing  thermodynamic  equation  of 
state  models  are  adequate  for  most  armor  penetration  calculations,  as 


are  existing  models  for  plasticity  (although  material  property  values 
related  to  plasticity  at  high  strain  rates  are  often  lacking).  The 
greatest  weakness  is  in  the  modeling  of  the  material  disintegration,* 
which  can  occur  by  several  mechanisms,  depending  on  the  material, 
temperature,  stress,  strain,  and  strain  rate.  Of  the  three  principle 
failure  modes  occurring  at  high  loading  rates — brittle  fracture,  ductile 
fracture,  and  shear  banding — the  last  two  have  been  most  commonly 
observed  in  connection  with  metallic  armor. 

Simple,  empirically  derived  failure  criteria  may  provide  some 
predictive  capability  for  penetration  scenarios  within  a  very  specific 
range  of  geometries  and  impact  velocities.  However,  a  generalized 
failure  model  applicable  to  a  wide  range  of  scenarios  requires  a 
detailed  qualitative  and  quantitative  understanding  of  the  material 
disintegration  processes  involved,  specifically  the  nucleation,  growth, 
and  coalescence  of  brittle  cracks,  ductile  voids,  and  shear  bands. 

Considerable  work  has  been  performed  in  studying  and  modeling  some 
of  these  material  disintegration  models.  SRI  International  has  devel¬ 
oped  the  nucleation  and  growth  (NAG)  Micro-Statistical  Fracture 

Mechanics  (MSFM)  models  for  brittle  and  ductile  dynamic  tensile  failure, 
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which  have  proven  quite  useful  as  predictive  models.  Before  this 

program,  some  work  had  been  done  in  studying  and  formulating  a  prelimi¬ 
nary  model  for  failure  by  the  nucleation,  growth,  and  coalescence  of 
11-14 

shear  bands.  However,  no  method  existed  that  would  incorporate 

these  models  into  a  computer  code  that  would  simulate  the  penetration  of 
armor  by  projectile  impact  and  predict  the  resulting  fragmentation 
environment. 


The  objective  of  this  program  was  to  develop  a  phenomenologically 
sound  material  disintegration  model  for  computationally  simulating  armor 
penetration  by  projectile  impact.  This  model  would  provide  the  capa¬ 
bility  for  predicting  the  downrange  fragment  size  and  velocity  distri¬ 
butions  and  would  eventually  be  used  by  the  DoO  for  armor  and  penetrator 
design. 


Our  approach  for  this  program  is  described  below. 

Baseline  materials — depleted  uranium  (DU)-3/4  wt%  titanium  alloy 
and  rolled  homogeneous  armor  (RHA) — were  procured  and  metallurgically 
analyzed.  Penetration  phenomenology  experiments  were  performed,  using 
various  impact  geometries,  velocities,  and  obliquities,  and  the  recov¬ 
ered  projectile  and  target  specimens  were  sectioned  and  examined  raetal- 
lographically  to  reveal  the  key  microstructural  damage  mechanisms  that 
govern  the  failure  and  disintegration  of  both  the  penetrator  and  the 
armor.  A  previously  acquired  data  base  was  available  for  4340  steel. 

Characterization  experiments  were  performed  that  exercised  a 
particular  failure  mechanism  in  a  geometry  that  is  relatively  simple  to 
simulate  computationally.  These  experiments  provided  a  basis  for 
developing  or  refining  a  failure  model  and  generating  material  failure 
kinetics  parameters.  The  contained  fragmenting  cylinder  (CFC)  test,  for 
example,  was  used  to  study  shear  band  nucleation  and  growth. 

Computational  models  were  developed  to  describe  the  observed 
failure  mechanisms.  The  models  were  based  upon  the  NAG  MSFM  models, 
which  attempt  to  qualitatively  and  quantitatively  simulate  the  actual 
micromechanical  failure  processes  involved. 

Material  failure  kinetics  parameters  were  determined  by  computa¬ 
tionally  simulating  the  material  characterization  experiments,  using  the 
failure  model  to  describe  the  damage.  The  parameters  were  then  used  to 
predict  damage  in  scenarios  with  different  geometries  or  loading 
conditions. 

A  computer  code  (C-HEMP)  was  assembled,  containing  the  features 
necessary  to  simulate  armor  penetration  by  projectile  impact.  Although 
this  code  includes  many  of  the  same  features  found  in  various  HEMP 
codes,  it  also  contains  new  logic  to  handle  the  complex  phenomenology 
observed  and  to  incorporate  the  material  failure  models  developed  above. 
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Finally,  computer  simulations  were  performed  to  assess  the  ability 
of  the  computational  method  Co  correctly  predict  the  main  features  of 
armor  penetration. 

These  six  parts  to  our  approach  for  this  program  were  not  performed 
sequentially,  but  proceeded  simultaneously  with  much  interaction  and 
iteration  between  them. 


Summary  of  Results  Obtained  During  First  Two  Years  of  Program 

Details  of  the  results  obtained  during  the  first  two  years  of  the 
program  are  available  in  the  first  and  second  annual  reports. 

Those  results  are  summarized  below. 

The  two  baseline  materials,  RHA  and  a  Dll  alloy,  were  obtained,  and 
RHA  was  examined  metallurgically .  Quasi-static  tensile  tests  were 
performed  on  RHA  specimens  that  were  cut  along  three  orthogonal 
directions.  The  resulting  stress-strain  flow  curves  revealed  a 
significant  anisotropy  in  ductility. 

We  studied  the  phenomenology  of  high-velocity  penetration  of 
moderate-to-thick  armor  plate  by  long-rod  projectiles  using  both  new  and 
existing  quarter-  and  full-scale  ballistics  data  for  DU  (and  4340  steel) 
projectiles  impacting  RHA  targets  at  velocities  both  near  and  well  above 
the  ballistic  limit,  and  at  obliquities  from  zero  to  seventy  degrees. 

We  drew  the  following  conclusions  from  these  studies: 

(1)  Shear  banding  is  the  dominant  failure  mechanism  in  both 
the  penetrator  and  the  target,  and  the  coalescence  of 
shear  bands  form  the  bulk  of  the  downrange  and  uprange 
fragments. 

(2)  Shear-band-induced  penetrator  backflow  and  shear  failure 
along  planes  of  weakness  in  an  anisotropic  target  are 
significant  mechanisms  of  target  erosion. 

(3)  Tensile  failure  may  be  a  factor,  but  only  in  a  rela¬ 
tively  few  instances,  such  as  at  the  rear  of  the  armor 
plate  in  the  final  stages  of  penetration  for  an  impact 
velocity  just  above  the  ballistic  limit. 
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(4)  Projectiles  change  direction  within  the  target  (oblique 
impacts  only)  seeking  the  path  of  least  resistance. 

(5)  Most  of  the  failure  mechanisms  in  a  highly  oblique 
impact  will  also  be  found  in  a  normal  impact  of  the  same 
materials  at  the  same  velocity. 

From  these  phemoraenological  conclusions,  we  selected  the  failure 
mechanisms  to  be  studied  in  detail  and  the  features  that  would  be  needed 
for  a  predictive  computer  code.  We  decided  that  shear  band  kinetics 
parameters  are  crucial  for  both  the  penetrator  and  target,  while  tensile 
failure  parameters  are  relatively  unimportant,  particularly  for  the 
penetrator.  Furthermore,  the  computer  code  would  need  features  that 
could  handle  the  extreme  shear  deformations  (a  rezoner),  the  rotation 
of  a  projectile  in  an  oblique  impact  (self -directed  or  wandering  slide 
line),  and  anisotropic  shear  failure.  Finally,  simulations  of  normal 
impacts,  which  are  easier  to  perform  than  those  of  oblique  impacts 
(because  of  the  axial  symmetry),  are  valid  tests  of  the  computational 
model's  ability  to  predict  f ragamentation  distributions. 

To  measure  the  kinetics  of  the  shear  banding  process,  we  performed 
(CFC)  experiments  with  RHA  specimens  cut  along  two  orthogonal  directions. 
The  experiments  yielded  a  wide  range  of  shear  banding  damage  and  confirmed 
the  expected  significant  anisotropy  in  shear  banding.  At  the  same  time, 
computational  model  for  shear  banding  in  an  isotropic  material  (SHEAR3) 
was  developed  and  calibrated  with  respect  to  previously  obtained  CFC 
data  on  4340  steel,  and  the  shear  band  kinetics  parameters  for  4340 
steel  (Rc40)  were  obtained. 

A  quarter-scale  armor  penetration  test  was  performed  at  BRL  to 
provide  a  back-of-the-armor  fragmentation  data  base  using  materials  (RHA 
and  4340  steel,  Rc40)  for  which  data  on  shear  band  kinetics  existed. 
Preliminary  simulations  of  this  test  were  made  using  the  SHEAR3  model 
(but  with  only  rough  estimates  of  the  RHA  shear  band  kinetics  param¬ 
eters)  and  an  existing  two-dimensional  code,  TROTT.  Although  the 
simulations  would  not  run  to  completion  (because  TROTT  could  not  handle 
the  large  shear  deformations),  the  calculated  size  distributions  of 
fragments  resulting  from  shear  banding  damage  were  in  fair  agreement 


with  experimental  observations,  indicating  that  a  phenomenologically 
based  failure  model,  which  describes  fragmentation  as  the  nucleation, 
growth,  and  coalescence  of  shear  bands,  has  a  good  potential  for  pre¬ 
dictive  capabilities. 

A  search  was  made  for  a  two-dimensional  wave  propagation  code  that 
had  all  or  most  of  the  computational  features  that  had  been  seen,  from 
the  phenomenology  study  described  above,  to  be  important  in  simulating 
armor  penetration  by  projectile  impact.  These  computational  features 
included  a  finite-element  cell  numbering  scheme  (to  allow  for  material 
ejection  and  for  the  inclusion  of  slide  lines  in  arbitrary  directions 
after  a  calculation  has  begun),  an  advanced  slide  line  capability 
(including  intersecting  and  wandering  slide  lines),  a  rezoning  feature 
that  can  account  for  extreme  cell  distortions,  and  the  ability  to  handle 
the  complex  failure  models.  When  this  search  turned  up  no  viable 
candidates,  work  was  begun  to  formulate,  assemble,  and  test  a  code 
Incorporating  many  of  these  features,  with  the  potential  for  adding  the 
remaining  features  at  a  later  time.  This  code  was  called  C-HEMP  (short 
for  Composite  HEMP),  because  it  was  largely  a  composite  of  several  HEMP- 
type  codes. 


Outline  of  Recent  Progress  and  Report  Organization 


In  this  section  we  will  briefly  outline  the  progress  made  in  this 
program  since  the  second  annual  report  and  describe  the  organization  of 
this  three-volume  report. 


The  experimental  work  planned  for  this  program  has  been  completed. 

A  metallurgical  analysis  was  performed  on  the  baseline  DU  alloy.  Quasi¬ 
static  tensile  tests  yielded  a  stress-strain  flow  curve,  and  the  DU 
specimens  recovered  from  these  tests  and  a  CFC  test  were  subjected  to 
metallographic  and  f ractographic  examination.  These  results  are  dis¬ 
cussed  in  Section  II  of  Volume  I. 


Quantitative  analysis  was  completed  on  the  RHA  CFC  experiments 
performed  during  the  first  and  second  years  of  the  program,  and  shear 
band  size  distributions  were  obtained.  A  series  of  DU  CFC  experiments 
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was  performed  and  analyzed,  and  the  shear  band  size  distributions  were 
obtained.  The  results  are  shown  in  Section  III  of  Volume  I,  along  with 
a  brief  review  of  the  CFC  technique. 

Refinements  were  made  to  the  SHEAR3  computational  model  for  shear 
banding  damage  in  isotropic  materials.  These  included  a  more  explicit 
strain  rate  dependence  for  shear  band  nucleation.  Complete  documenta¬ 
tion  for  this  model  is  in  Volume  II  of  this  report.  Work  was  begun  to 
develop  a  significantly  more  complex  model  for  shear  banding  in  aniso¬ 
tropic  materials  (SHEAR4),  but  limitations  on  the  scope  of  this  program 
did  not  allow  this  work  to  be  completed.  A  discussion  of  some  of  the 
features  of  SHEAR4  is  included  in  Section  VI  of  Volume  I. 

The  modified  version  of  SHEAR3  was  used  (in  conjunction  with 
C-HEMP)  to  simulate  the  RHA  CFC  experiments,  to  obtain  a  preliminary  set 
of  shear  band  kinetics  parameters  for  RHA,  and  to  obtain  a  revised  set 
of  parameters  for  4340  steel  (Rc40).  The  results  are  shown  in  Section 
IV  of  Volume  I. 

The  development  and  testing  of  the  C-HEMP  code  is  nearly  complete. 
Incorporated  features  that  have  proven  successful  in  a  several  different 
applications  include  a  free-field  input  format,  a  finite-element  cell 
numbering  scheme,  a  rezoner,  a  reasonably  complex  slide  line  routine 
(including  slip,  with  or  without  a  frictional  component,  along  user- 
located  single  or  intersecting  slide  lines  that  are  either  already 
"unzipped"  or  "unzip"  as  the  calculation  proceeds),  and  the  inclusion  of 
the  SHEAR3  and  elastic-plastic  equation-of-state  subroutines.  Complete 
documentation  of  C-HEMP  will  be  available  later. 

A  series  of  computational  simulations  of  armor  penetration  was 
performed  with  C-HEMP,  to  test  the  various  features  of  the  new  code. 


These  simulations  of  the  normal  impact  of  a  4340  steel  projectile  into 
an  RHA  armor  plate  demonstrate,  at  least  in  a  preliminary  way,  the 
ability  of  C-HEMP  and  SHEAR3  to  realistically  model  armor  penetration 
both  near  and  well  above  the  ballistic  limit.  This  work  appears  in 
Section  V  of  Volume  I. 

Finally,  we  have  written  a  detailed  set  of  recommendations  for 
future  work  in  the  area  of  computational  modeling  of  armor  penetra¬ 
tion.  This  appears  in  Section  VI  of  Volume  I. 


II  METALLURGICAL  ANALYSIS  OF  DU 
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Computational  modeling  of  the  large  deformations  and  material 
failure  that  accompany  the  penetration  of  an  armor  plate  by  a  projectile 
requires  knowing  the  constitutive  behavior  of  the  materials  involved  and 
identifying  the  microstructural  mechanisms  leading  to  fracture.  Much  of 
this  information  can  be  provided  by  quasi-static  tensile  tests  and 
metallographic  and  f ractographic  observations  of  the  tensile  specimens. 

Metallurgical  studies  of  the  baseline  armor  material,  RHA,  conducted 
during  the  first  year  of  the  program,  have  been  previously  reported.^ 
Similar  studies  were  performed  during  the  final  year  of  the  program  on 
the  baseline  penetrator  material,  a  DU-3/4Ti  alloy.  In  this  section,  we 
will  discuss  the  results  of  these  studies  and  a  further  metallurgical 
study  of  a  DU  specimen  recovered  from  a  CFC  experiment. 


Duasistatic  Tensile  Tests 


Smooth  round-bar  tensile  specimens  were  machined  by  the  Rocky  Flats 
Plant  of  Rockwell  International  in  Golden,  CO,  from  extruded  and  heat- 
treated  rods  of  a  single  heat  of  DU-3/4Ti  alloy.  Three  sets  of 


specimens  were  prepared  in  accordance  with  specifications  provided  by 
SRI.  The  first  set  of  specimens  had  a  2.5-cm  gage  length,  a  0.625-cm 
gage  diameter,  and  was  machined  with  the  axes  parallel  to  the  extrusion 
direction  (longitudinal  specimens).  The  second  set  had  a  1.27-cm  gage 
length,  a  0.32-cm  diameter,  and  was  also  machined  with  the  axes  perpen¬ 
dicular  to  the  extrusion  direction  and  intersecting  the  rod  centerline 
(transverse  specimens).  The  third  set  had  a  0.51-cm  gage  length,  a 
0.32-cm  gage  diameter,  and  was  also  machined  in  the  transverse  direc¬ 
tion,  but  not  intersecting  the  rod  centerline  (microtensile  specimens). 
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This  last  set  was  for  the  purpose  of  assessing  the  affect  of  any 

* 

centerline  porosity  on  the  tensile  properties  of  the  extruded  rods. 

Two  longitudinal  and  two  transverse  tensile  specimens  were  pulled 
to  fracture,  using  an  Instron  testing  machine  at  a  crosshead  speed  of 
8  x  10~4  cm/s.  The  results  of  these  tests  are  summarized  in  Table  l 
and  the  corresponding  true  stress-true  strain  curves  are  shown  in 
Figure  1. 

The  results  show  that  in  the  the  DU  alloy  is  somewhat  weaker  in 
the  transverse  direction  than  in  the  longitudinal  direction.  For  trans¬ 
verse  directions  the  yield  stress  is  about  7.6  kbar  (110  ksi)  and  true 
fracture  stress  is  about  15  kbar  (218  ksi),  as  compared  to  9.0  kbar 
(130  ksi)  and  16.2  kbar  (235  ksi),  respectively,  for  the  longitudinal 
specimens.  On  the  other  hand,  the  ductility  is  about  the  same  for  both 
specimen  orientations,  considering  the  scatter  in  the  elongation  and 
fracture  strain  data.  This  observation  is  consistent  with  the  litera¬ 
ture  on  the  dependence  of  fracture  toughness  on  orientation  for  extruded 
DU-3/4T1.1  From  these  results,  we  concluded  that  any  centerline  poro¬ 
sity  that  might  exist  is  not  important  to  the  mechanical  properties  of 
these  rods,  and  therefore  testing  of  the  microtensile  specimens  was 
unnecessary. 

The  post-yield  data  in  Figure  1  may  be  described  by  a  power  law 
hardening  equation  of  the  form: 


where  k  =  24.9  kbar  (236  ksi)  and  n  «*  0.20  for  the  longitudinal  data, 
and  k  =  23.5  kbar  (234  ksi)  and  n  =*  0.22  for  the  transverse  data. 


Centerline  microporosity  has  been  a  chronic  problem  in  the  extrusion 
of  uranium  alloys.  Although  the  pores  are  often  quite  small  (smaller 
than  the  detection  limit  of  x-rays),  they  have  the  potential  for  influ¬ 
encing  material  properties  of  extruded  products.  The  cause  of  such 
porosity  is  not  well  established. 
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Table  1 


SUMMARY 

OF  QUASI-STATIC 

TENSILE 

TEST 

RESULTS  ON  DU-3/4 

Ti  ALL l)Y 

Specimen 

Yield 

(0.2% 

stress 

Offset) 

True  Stress 
at  Fracture 

True  Strain 

Elongation  to 

Number 

kbar 

ksi 

kbar 

ksi 

at  Fracture 

Failure  {l) 

Longitudinal 

06R31 

8.97 

130 

16.5 

239 

0.155 

17 

06R32 

8.97 

130 

15.9 

231 

0.115 

12 

Transverse 

06MTI 

7.73 

112 

15.4 

223 

0. 148 

16 

06MT2 


7.45  108 


14.7  213 
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FIGURE  1  MEASURED  STRESS-STRAIN  CURVE  OF  DU-3/4  Ti  ALLOY 
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TRUE  STRESS,  a  (kbar) 


These  minor  differences  in  work-hardening  parameters  are  most  likely  due 
to  microstructural  texturing  effects  during  the  extrusion  process. 


Metallography 

The  broken  halves  of  one  longitudinal  and  one  transverse  DU  tensile 
specimen  were  sectioned  along  their  axis,  ground,  and  electropolished 
for  metallographic  examination.  Figure  2  is  a  micrograph  from  the 
longitudinal  specimen.  The  needle-like  character  of  the  raicrostructure 
is  due  to  the  martensitic  (diffusionless)  transformation  of  the  high 
temperature  q-phase  (BCC)  into  the  a-phase  (orthorhombic)  upon  quenching. 
The  as-quenched  structure  is  not  equilibrium  a-phase,  but  is  rather  a 
slightly  distorted  orthorhombic  structure  termed  a'.  The  q  grain 
structure  generated  during  the  solutionizing  treatment  is  still  faintly 
visible  in  Figure  2.  The  q  grain  size  is  estimated  to  have  been  on  the 
order  of  0.1  to  0.5  mm. 

The  martensitic  transformation  of  q  uranium  into  the  metastable 
structure  did  not  strengthen  the  material  (as  it  would  in  steel);  it  is 
simply  a  mechanism  to  achieve  a  supersaturated  uranium-titanium  solid 
solution.  The  material  was  strengthened  by  a  subsequent  aging  treatment 
at  360°C,  as  will  be  discussed  below.  It  is  evident  from  Figure  2  that 
this  aging  treatment  did  not  substantially  alter  the  acicular  morphology 
of  the  a'  grains,  suggesting  that  only  a  small  amount  of  the  raetas table 
a'  converted  to  equilibrium  a  during  aging.  Although  not  visible  in 
Figure  2,  patches  of  equilibrium  a  were  occasionally  observed  along  the 
a  grain  boundaries. 

Note  also  from  Figure  2  that  the  alloy  contains  many  large  (approx¬ 
imately  10  urn  across)  precipitates.  Energy-dispersive  x-ray  analysis 
indicates  that  these  precipitates  are  Ti  rich,  and  probably  consist  of 
the  intermetallic  compound  U2H.  In  longitudinal  cross  sections — that 
is,  in  planes  containing  the  rod  axis — these  precipitates  are  in  ordered 
long  stringers  parallel  to  the  extrusion  direction.  In  transverse  cross 
sections,  the  precipitate  distribution  is  random.  Therefore,  these 
precipitates  form  either  during  cooling  from  the  homogenizing  temperature 
prior  to  extrusion  or  during  the  extrusion  process  itself. 


In  addicion  Co  these  large  unintentional  precipitates,  there  is  a 
uniform  distribution  of  very  fine  («  1  pm  across)  l^Ti  particles,  not 
resolvable  by  optical  microscopy.^  These  particles  are  produced  by 
aging  the  solutionized  and  quenched  extruded  rod  at  360°C  for  several 
hours,  and  they  inhibit  the  motion  of  dislocations  through  the  lattice. 
Therefore,  a  higher  stress  is  needed  to  move  the  dislocations  past  these 
precipitates,  increasing  the  yield  strength  of  the  material. 

Examination  of  the  cross  sections  of  the  broken  tensile  specimens 
did  not  reveal  any  microcracking  or  void  growth  below  the  fracture 
plane.  Furthermore,  all  the  specimens  failed  before  any  significant 
necking  occurred.  However,  the  second  longitudinal  specimen  (06R32  in 
Figure  1)  did  display  an  interesting  feature:  a  secondary  crack  formed 
perpendicular  to  the  axis  and  about  0.5  cm  below  the  main  fracture 
plane.  This  crack  propagated  through  95%  of  the  cross  section  before 
arresting. 

Fractography 

The  fracture  surfaces  of  a  longitudinal  and  transverse  DU  tensile 

specimen  were  studied  with  scanning  electron  microscopy.  In  both  cases, 

the  surface  was  a  mixture  of  transgranular  and  dimple  fracture  as  has 
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been  observed  in  the  literature.  The  surface  topology  of  a  repre¬ 
sentative  longitudinal  specimen  is  illustrated  in  Figures  3  and  4. 

Figure  3  shows  large  steps  in  the  fracture  surface  running  from  bottom 
left  to  top  right  across  the  fractograph.  The  darker  featherlike 
features  are  the  martensitic  needles  seen  previously.  Figure  4  shows  a 
detail  of  one  of  the  steps  in  the  fracture  surface,  with  secondary 
cracks  running  parallel  to  the  cleavage  plane.  Note  also  the  network  of 
dimples  surrounding  the  transgranular  fracture. 

Several  large  globular  U2Ti  precipitates  are  also  visible  in 
Figure  3.  While  these  precipitates  are  too  large  to  affect  the  yield 
strength,  they  definitely  contribute  to  the  ultimate  tensile  fracture 
process  and  lower  the  overall  ductility. 
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The  only  difference  between  the  fracture  surface  of  the  longitudinal 
and  transverse  specimens  is  the  already  mentioned  difference  in  the  large 
precipitate  distribution.  On  the  fracture  surface  (which  is  perpendicular 
to  the  cross  sections  used  for  metallographic  observations),  the  preci¬ 
pitate  distribution  is  random  for  the  longitudinal  specimens,  whereas 
long  ordered  stringers  are  seen  in  the  transverse  specimens.  A  typical 
large  U2H  precipitate  is  shown  in  Figure  5.  It  has  a  characteristic 
cubic  shape  and  is  about  10  pm  in  across.  The  matrix-l^Ti  interface  is 
apparently  stronger  than  the  l^Ti,  and  thus  in  most  cases  the  inclusions 


fracture  rather  than  debond. 


Metallurgical  Analysis  of  DU  CFC  Specimen 


As  will  be  discussed  in  detail  in  Section  III,  we  performed  a 
series  of  DU  CFC  experiments  to  study  shear  banding  kinetics  in  DU.  T 
CFC  specimens  were  cylindrical  tubes  machined  from  extruded  rod  from 
the  same  heat  used  for  the  tensile  specimens,  and  were  thermoraechani- 
cally  processed  in  a  similar  manner.  We  metallurgically  examined  a  DU 
specimen  recovered  from  one  of  these  experiments  to  determine  if  the 
explosively  driven  deformation  and  heating  that  accompany  a  CFC  experi 
ment  caused  any  significant  microstructural  changes  in  the  material. 


The  microstructure  of  the  recovered  DU  CFC  fragment,  as  shown  in 
Figure  6(a),  is  primarily  acicular  ct',  the  martensitic  form  of  the 
equilibrium  orthorhombic  (a-U)  structure.  This  martensitic  form, 


produced  by  rapidly  quenching  U-Ti  alloys  from  the  BCC  ( y)  phases,  is 


also  orthorhombic,  but  with  a  shortened  b  axis.  The  prior  y  grain 


boundaries,  clearly  visible  in  Figure  6(a),  are  the  result  of  the  pre¬ 


quench  soaking  treatment  at  temperatures  higher  than  750°C. During 
quenching,  the  a'  platelets  initiate  from  the  y  grain  boundaries  and 
propagate  rapidly  across  the  grains,  thereby  marking  the  locations  of 
these  grain  boundaries.  This  is  essentially  the  same  microstructure  as 
observed  for  the  tensile  specimens. 


V  N."  /aC 


FIGURE  5  U2Ti  INCLUSION  AND  DUCTILE  DIMPLES 
ON  FRACTURE  SURFACE  OF  DU-3/4  Ti 
ALLOY 
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Note  also  In  Figure  6(a)  Che  white  material  found  predominantly 
along  the  prior  gamma  grain  boundaries.  This  material  is  the  equi- 
librium  a-U/U^Ti  structure  and  is  a  result  of  overaging.  During 
prolonged  aging  treatments,  the  equilibrium  structure  nucleates  along 

the  prior  boundaries  and  then  grows  into  the  matrix.  If  held  long 

★ 

enough  at  temperature,  Che  entire  microstructure  would  revert  to  the 
equilibrium  structure  of  a-U/U^Ti.  As  shown  in  Figure  b(b),  this  equi¬ 
librium  structure  initially  e>  ibits  a  lamellar  morphology  with  alter¬ 
nating  layers  of  oc-U  and  l^Ti.  If  held  for  very  long  times  at  temper¬ 
ature,  the  L^Ti  platelets  will  eventually  become  spheroid  to  minimize 
the  free  energy  of  the  system. 

Since  both  the  aclcular  a'  structure  and  the  equilibrium  a-U/U  Ti 
structure  were  found  in  the  quasi-static  tensile  specimens,  we  concluded 
that  the  shock  loading  and  subsequent  heating  did  not  substantially 
alter  the  microstructure.  To  confirm  this  we  performed  a  series  of 
hardness  tests  across  the  thickness  of  the  recovered  fragment  and 
measured  a  hardness  of  395  ±  16  VHN.  This  differed  negligibly  from  the 
388  ±  26  VHN  hardness  measured  in  the  as-received  specimens.  It  there¬ 
fore  appears  that  the  CFC  experiment  does  not  cause  any  significant 
microstructural  changes  that  may  influence  the  initiation  or  propagation 
of  adiabatic  shear  bands. 

Summary  of  Results 

The  DU-3/4T1  alloy  from  extruded  rod,  which  was  tested  for  this 
program,  has  the  following  tensile  properties: 

yield  stress  —  7. 4-9.0  kbar  (108-130  ksi) 
true  fracture  stress  —  14.7-16.5  kbar  (213-239  ksi) 
true  fracture  strain  —  12-16% 


Unlike  the  martensitic  a — *-Ot '  transformation,  the  a' — *a  transfor¬ 
mation  is  diffusion  controlled. 
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Relatively  little  anisotropy  in  tensile  properties  was  observed  between 
longitudinal  and  transverse  directions. 

The  DU  alloy  has  an  acicular  a'  microstructure  with  a  nonuniforra 
distribution  of  large  U2H  precipitates.  Fracture  occurs  by  a  mixture 
of  cleavage  and  dimple  rupture,  with  cleavage  as  the  dominant  mode. 

No  siginificant  microstructural  changes  were  found  in  rapidly 
deformed  specimens,  before  initiation  of  shear  banding. 


Ill  MATERIAL  CHARACTERIZATION  EXPERIMENTS 


During  this  program,  we  performed  two  series  of  CFC  experiments, 
one  with  RHA  and  one  with  DU.  The  purpose  of  these  experiments  was  to 
measure  the  shear  band  size  and  density  distributions  in  specimens  that 
underwent  a  readily  characterized  strain  history,  and  determine  from 
this  data  the  shear  banding  kinetics  parameters. 

We  have  previously  reported^-^  details  of  the  experimental  tech¬ 
nique  and  the  analytical  method  for  obtaining  shear  band  kinetics  param¬ 
eters  from  the  experimental  data  (as  was  performed  last  year  for  4340 
steel,  Rc40).  In  this  section  we  will,  after  a  brief  review  of  the  CFC 
technique,  present  the  results  of  the  CFC  experiments  on  RHA  and  DU. 

The  CFC  technique  consists  of  applying  internal  explosive  loads  to 
thick-walled  cylinders  of  the  specimen  material  so  that  they  expand 
rapidly  until  stopped  by  massive  concentric  containment  cylinders.  The 
setup  is  shown  schematically  in  Figure  7.  The  initial  deformation  rate 
and  the  final  deformation  achieved  are  controlled  by  the  explosive  den¬ 
sity  and  the  thicknesses  of  the  various  containment  cylinders.  Mode  L 
shear  bands  (see  Figure  8)  originating  at  the  inner  surface  of  the 
recovered  specimen  tube  are  measured  and  counted,  and  the  resulting 
shear  band  size  distributions  as  a  function  of  axial  position  along  the 
tube  are  correlated  with  those  obtained  from  computer  simulations  to 
determine  the  shear  banding  kinetics  parameters. 


CFC  Experiments  with  RHA 


Nine  RHA  CFC  experiments  were  performed  during  the  second  year  of 
the  program.  The  experimental  configuration  and  explosive  parameters 
for  these  experiments  are  shown  in  Figure  9  and  in  Table  2,  along  with 
the  qualitative  damage  results. 


Because  quasi-static  tensile  tests,  performed  during  the  first  year 
of  the  program,  had  revealed  a  significant  anisotropy  in  tensile  failure 
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FIGURE  7 


SCHEMATIC  OF  CFC  EXPERIMENTS  FOR  STUDYING 
SHEAR  BAND  KINETICS 
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8  —  Length  in  Direction 
Perpendicular  to 
Slip  Motion 

d  —  Depth  Along  45& 
Slip  Plane 

B  -  Shear  Displacement 
Along  Slip  Plane 

Z  —  Axial  Direction 

r  —  Radial  Direction 

8  —  Circumferential 
Direction 


ORIENTATION  MODE  3 
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FIGURE  8  GEOMETRY  AND  NOMENCLATURE  FOR  SHEAR  BANDS 
IN  FRAGMENTING  CYLINDER  EXPERIMENTS 
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FIGURE  9  CROSS  SECTION  OF  RHA  CFC  EXPERIMENTS 


Table  2 


PARAMETERS  AND  QUALITATIVE  RESULTS  FOR  RHA  CFC  EXPERIMENTS 


Specimen 

Lucite 

Buffer 

Thickness 

(cm) 

PETN 

Density 

(g/cm3) 

Calculated 

Detonation 

Pressure 

(kbar) 

Qualitative  Shear 

Banding  Results 

L-4* 

0.873 

1.2 

1 22 

Substantial  Mode  1  banding  con¬ 
centrated  in  rolling  direction 

L-2 

0.73 

1.0 

80 

Same  as  L-4 

t-it 

0.873 

1.2 

122 

A  few  incipient  Mode  1  bands 
(some  tensile  cracking  on  outer 
surface  in  midaxial  region). 

T-2 

1.600 

1.2 

122 

Few  incipient  Mode  1  bands. 
Moderate  Mode  2  banding  in  mid- 
axial  region. 

T-3 

0.873 

1.35 

162 

Incipient  Mode  1  banding  over 
most  of  inner  surface  (some 
tensile  cracking  as  in  T-l 
above) . 

T-5 

1.600 

1.35 

162 

Extensive  incipient  Mode  1 
banding.  Few  moderate  Mode  2 
bands  in  midaxial  region. 

T-7 

1.600 

1.47  to 
»  1.77 

202  to 
-  324 

Substantial  Mode  1  banding. 
Substantial  Mode  2  banding  in 
narrow  midaxial  region. 

T-9 

1.27 

1.54 

228 

Substantial  Mode  1  banding. 

Mode  2  bands  in  midaxial  region 
cuts  specimen  in  two. 

T- 10 


1.27 


1.6 


250 


Same  as  T-9,  with  some 
f  ragmentat ion 


properties  in  RHA,^  we  expected  a  similar  anisotropy  to  exist  with 
respect  to  shear  banding.  Early  experiments  with  RHA  specimen  tubes 
whose  axes  were  in  the  long  transverse  direction  of  the  RHA  slab  (the 
L-specimens,  as  shown  in  Figure  10)  confirmed  this  expectation.  The 
preponderance  of  shear  bands  in  planes  parallel  to  the  rolling  plane  in 
specimens  L-2  and  L-4  (as  seen  in  Figures  11  and  12),  for  example, 
clearly  showed  that  shear  banding  will  occur  at  lower  strains  in  the 
rolling  plane  than  in  planes  perpendicular  to  the  rolling  planes. 

Once  it  became  clear  that  by  using  L-specimens  we  could  not  obtain 

a  uniform  distribution  of  shear  bands  around  the  circumference  of  the 

specimen,  we  switched  to  T-specimens  (whose  axes  were  in  the  through- 

the-thickness  direction)  for  which  all  the  Mode  1  bands  would  lie  in 

planes  perpendicular  to  the  rolling  plane.  The  first  four  T-specimen 

experiments  (T-l,  T-2,  T-3,  and  T-5),  using  PETN  explosive  densities 
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similar  to  those  in  previous  experiments,  *  yielded  only  incipient 
Mode  1  shear  banding.  However,  the  final  three  experiments  (T-7,  T-9 , 
and  T-10),  using  much  higher  explosive  densities,  yielded  significant 
quantities  of  both  incipient  and  large  shear  bands,  as  well  as  some 
fragmentation.  Photographs  of  these  specimens  are  shown  in  Figures  13 
and  14. 

Shear  band  size  distributions  were  obtained  from  the  three  RHA 
specimens  that  exhibited  significant  amounts  of  Mode  1  shear  banding  in 
planes  perpendicular  to  the  rolling  plane.  (No  attempt  was  made  to 
obtain  quantitative  damage  estimates  for  the  few  Mode  2  bands  found  in 
these  T-specimens,  or  for  the  Mode  1  bands  in  the  L-specimens  located  in 
the  rolling  plane.  These  few  bands  localized  in  narrow  regions  would 
have  yielded  poor  statistical  results. )  The  cumulative  size  distri¬ 
butions  in  the  various  axial  zones  in  shots  T-7,  T-9,  and  T-10  are  shown 
in  Figures  15,  16,  and  17,  respectively.  The  areal  density  of  the  shear 
bands  intersecting  the  inner  surface  of  the  specimen  tube  and  having  a 
length  (as  defined  in  Figure  8)  greater  than  L,  is  plotted  against  L. 

As  expected,  for  each  shot  the  zones  at  or  near  the  ends  of  the  specimen 
exhibited  lower  shear  band  densities  than  those  in  the  middle. 
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FIGURE  10  CFC  SPECIMENS  MACHINED  FROM  RHA  SLAB,  SHOWING  MODE  1  SHEAR  BANDS 
IN  PLANES  PERPENDICULAR  TO  ROLLING  PLANE 


FIGURE  12  POLISHED  AND  ETCHED  SECTION  THROUGH  RHA  SPECIMEN  TUBE  L-4 
SHOWING  PREPONDERANCE  OF  SHEAR  BANDING  DAMAGE  IN  ROLLING 
DIRECTION 

/Plane  s  cut  at  a  distance  of  2-3  8"  from  bottom  of  specimen  tube.) 
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SHEAR  BAND  LENGTH,  L  (cm) 

MA-7893-1  51 

FIGURE  15  CUMULATIVE  SIZE  DISTRIBUTIONS  OF  SHEAR  BANDS  FROM  RHA  SHOT  T-7 
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SURFACE  DENSITY  OF  SHEAR  BANDS  WITH  LENGTH 


There  Is  a  significant  uncertainty  in  the  plotted  values  for  the 
shear  band  density  for  lengths  greater  than  about  0.7  cm  because  the 
small  number  of  these  large  bands  produce  poor  statistics.  There  Is 
also  a  significant  uncertainty  in  the  plotted  values  for  lengths  less 
than  0. L  cm,  because  it  was  difficult  to  see  and  accurately  count  the 
large  number  of  small  bands  (this  was  particularly  true  of  RHA  because 


the  low  carbon  content  minimizes  the  martensitic  transformations  that 


appear  as  white  etching  bands).  So  the  distributions  are  most  accurate 
in  the  intermediate  length  range. 


CFC  Experiments  with  DU 


A  series  of  three  CFC  experiments  were  performed  using  DU.  The 
experimental  configuration  and  explosive  parameters  for  these  experi¬ 
ments  are  shown  in  Figure  18  and  Table  3,  along  with  the  qualitative 
damage  results.  Note  that  the  buffer  thickness  was  kept  constant  for 
all  three  shots,  but  the  PETN  density  and  hence  the  detonation  pressure 


was  varied. 


Predictions  from  Critical  Strain  Theory 


We  expected  that  DU  would  exhibit  shear  bands  at  lower  plastic 


strains  than  those  in  the  other  materials  (RHA  or  4340  steel,  R^Q)  with 


which  we  have  performed  CFC  studies.  This  expectation  was  based  on 


calculations  using  our  previously  derived  formula1  for  critical  strains 


for  shear  band  initiation: 


2pEn  (n  +  1)  \n+l 
- 2 - 

n  ( 2n  +  1 )  / 
o  / 


where 


v  s  equivalent  plast  strain  at  which  shear  bands  mav  nucleate 
cr 


0  H  density 


E  =  specific  internal  energy  at  incipient  melt. 


and  n  are  defined  by  the  strain-hardening  formula 
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CROSS  SECTION  OE  DU  CEC  EXPERIMENTAL  APPARATUS 


Table  3 
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PARAMETERS  AND  QUALITATIVE  RESULTS  FOR  DU  CFC  EXPERIMENTS 


Shot 

PETN 

Density 

p  (g/cm3) 

Explosive 
Column 
Height, 
h  (cm) 

Aery  Lie 
Buffer 
Thickness 

(cm) 

Expected  Peak 
Detonation 

Pressure 
(kbar ) 

Qualitative  Damage 
Results 

10 

1.20 

7.62 

0.762 

120 

Complete  fragmenta¬ 
tion  by  Mode  1  shear 
banding 

11 

0.95 

7.62 

0.762 

74 

Substantial  Mode  l 
shear  banding, 
moderate  fragmenta¬ 
tion. 

12 

0.80 

6.67 

0.762 

45 

Substantial  Mode  1 
shear  banding,  minor 
fragmentation 

39 
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where 


r  =  shear  scress 


£ ?  =  equivalent  plastic  strain 
and  a  is  defined  by  the  thermal  softening  function 


f/M-  i  -  & 

E  E 

\  m  /  \  m 


where 


W  =  specific  energy  of  the  plastic  deformation. 

Table  4  gives  the  claculated  critical  strains  for  various  materials, 
including  DU,  along  with  the  material  parameters,  as  defined  above,  that 
were  used  in  the  calculation.  The  critical  strain  calculated  for  DU  is 
0.10,  as  compared  with  0.56  for  RHA  and  0.19  for  4340  steel,  Rc40. 
Previous  CFC  experiments  on  RHA  and  4340  steel  have  shown  this  formula, 
which  was  derived  from  plastic  instability  theory,  to  be  useful  in 
predicting  approximate  strains  for  shear  band  initiation. 


Qualitative  Experimental  Results 


Examinations  of  the  recovered  specimens  showed  that  our  predictions 
were  reasonable.  Figure  19  shows  photographs  of  the  two  tubes  that  did 
not  completely  fragment  (shots  11  and  12),  sliced  in  half  axially  to 
allow  easier  inspection.  Although  the  shear  bands  on  the  inner  surface 
are  somewhat  obscured  by  the  substantial  oxidation  or  scale,  it  appears 
Chat  for  shot  12,  Mode  1  shear  bands  appear  at  radial  expansions  of 
about  10%. 

As  seen  in  Figure  19(a),  a  transverse  cut  near  the  midaxial  plane  was 
made  in  one  DU  piece  recovered  from  shot  11.  The  ground  and  electro- 
polished  surface  (in  Figure  20)  shows  several  large  cracks  and  several 
smaller  bands  crossing  the  specimen  at  approximately  45°  to  the  radial 
direction.  Meta  1 lographic  examination  showed  that  both  the  cracks  and 
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Table  4 

CRITICAL  STRAINS  FOR  SHEAR  BAND  INITIATION 
IN  VARIOUS  MATERIALS,  AS  CALCULATED  FROM  EQUATION  (L), 
AND  MATERIAL  PARAMETERS  USED  IN  THE  CALCULATIONS 


Material 


e15  (dyne/cm^) 

cr 


(ergs/ 4) 


ote:  Numbers  in  parentheses  refer  to  the  sources  for  the  data,  which 

listed  in  the  references. 

Best  estimates  from  available  data  on  similar  materials. 

Data  obtained  in  this  program — see  Section  II. 
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Che  bands  are  associated  with  Che  adiabatic  shearing  process.  Higher 
magnification  views  of  two  of  the  shear  bands  (Figure  21)  reveals  a 
white,  featureless  zone  bounded  by  a  heavily  sheared  microstructure,  as 
evidenced  by  Che  curvature  of  the  martensite  platelets  adjacent  to  the 
band.  This  structure  is  similar  to  that  observed  previously  in  a  DU 
penetrator  fragment.^ 

Anisotropy 

The  CFC  experiments  are  designed  to  produce  large  quantities  of 
Mode  I  bands  because  the  Mode  1  bands  lie  in  planes  of  maximum  shear 
strain  (parallel  to  the  cylinder  axis  and  at  45°  from  the  radial  and 
circumferential  directions).  The  DU  CFC  specimens  were  machined  from 
extruded  rod  in  the  only  orientation  feasible,  that  is,  with  the 
specimen  axis  co-axial  with  the  rod.  The  Mode  L  shear  band  orientation 
produced  in  the  CFC  experiments  is,  unfortunately,  not  the  same  orien¬ 
tation  as  bands  primarily  responsible  for  penetrator  erosion  during 
armor  penetration.  Previous  phenomenological  studies  with  several 
different  penetrator  materials ^ ^ ^  have  shown  that  such  erosion  is 
caused  by  Mode  2  bands:  bands  parallel  to  the  circumferential  direction 
and  at  45“  from  the  axial  and  radial  directions. 

There  may  be  a  few  Mode  2  bands  in  the  recovered  DU  specimens.  For 
example,  the  shear  lip  at  the  extreme  right  edge  of  Figure  13(b)  is  a 
Mode  2  band  caused  by  edge  effects.  However,  our  experience  indicates 
that  in  a  CFC  experiment  only  Mode  L  bands  are  produced  in  large  enough 
numbers  and  over  a  wide  enough  region  to  yield  quantifiable  results. 
Therefore,  we  will  have  to  use  the  results  of  the  Mode  1  analysis  to 
predict  Mode  2  damage.  This  might  cause  difficulties  in  the  case  of 
significant  shear  banding  anisotropy. 

We  have  two  significant,  although  not  totally  conclusive,  pieces  of 
evidence  indicating  that  anisotropy  is  not  a  major  factor  in  this  UU 
alloy.  First,  there  was  negligible  Mode  2  shear  banding  observed  in  the 
recovered  CFC  specimens.  This  suggests  that  planes  in  the  Mode  2  orien¬ 
tation  in  DU  are  not  much  weaker  in  their  resistance  to  shear  banding 
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than  those  in  the  Mode  1  orientation.  In  contrast,  RHA,  which  was  found 
to  exhibit  just  such  a  shear  banding  anisotropy,  did  show  significant 
Mode  2  shear  banding  in  recovered  CFC  specimens  whose  axes  were 
oriented  transverse  to  the  rolling  plane  (see  Figure  14,  for  example). 

Second,  static  tensile  tests  performed  on  DU  (as  reported  in  the 
Section  II)  showed  only  a  negligible  anisotropy  between  the  axial  and 
transverse  directions  for  failure  strain.  We  believe  the  absence  of 
anisotropy  in  tensile  failure  is  a  good  predictor  of  its  absence  in 
shear  failure,  since  anisotropy  in  both  failure  modes  is  thought  to 
arise  from  the  same  microstructural  features. 

Quantitative  Damage  Analysis 


Of  the  three  DU  CFC  experiments  performed,  only  one — shot  12  [as 
shown  in  Figure  9(b)] — appeared  to  have  achieved  a  sufficiently  low 
level  of  shear  banding  to  allow  for  a  quantitative  damage  assessment. 
However,  the  inside  surfaces  of  the  recovered  tube  was  covered  by  a 
thick  black  coating,  undoubtedly  the  result  of  the  reaction  between  the 
pyrophoric  uranium  and  the  explosive  products.  Because  the  thick 
coating  obscured  much  of  the  detail  of  the  shear  banding  damage,  it  was 
necessary  to  remove  the  coating  before  proceeding  with  the  analysis. 

After  several  different  chemical  cleaning  techniques  (based  on 
established  procedures  for  removing  uranium  oxide)  proved  unsatisfac¬ 
tory,  we  tried  an  electrochemical  technique  recommended  by  personnel  at 
the  Lawrence  Livermore  Laboratory.  This  technique  used  a  solution  of 
50  txL  of  perchloric  acid,  200  mL  of  acetic  acid,  and  2  grams  of  citric 
acid,  cooled  by  a  liquid  nitrogen  bath.  A  current  density  of  100  to  150 
raA/cm1"  was  applied  to  the  specimen  for  a  few  minutes,  using  a  stainless 
steel  cathode  immersed  in  the  solution  to  complete  the  circuit.  The 
electrochemical  technique  removed  all  the  scale  except  for  a  thick  band 
that  was  subsequently  shown  to  be  a  carbon-based  material,  rattier  than 
uranium  oxide.  Fortunately  this  band  was  easily  removed  bv  brushing  mu 
scraping  the  surface. 
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The  cleaned  specimen  Is  shown  in  Figure  22,  along  with  the  axial 
zones  used  in  the  quantitative  shear  band  damage  analysis.  The  number 
of  shear  bands  of  various  lengths  were  counted  for  each  zone,  and  the 
resultant  cumulative  size  distributions  are  shown  in  Figure  23.  Note 
that  the  uncertainty  in  these  surface  density  curves  is  greater  than 
that  for  RHA  or  4340  steel.  The  condition  of  the  DU  surface,  even  after 
cleaning,  tended  to  obscure  bands  less  than  a  couple  of  mm  in  length  and 
to  make  it  difficult  to  determine  if  a  long  band  is  one  band,  or  two  or 
more  joined  together. 

The  quantitative  shear  band  results,  along  with  the  measurements  of 
final  wall  thickness  and  diameter  as  a  function  of  axial  position,  are 
therefore  now  available  for  current  and  future  comparison  with  computer 
simulations  to  determine  the  shear  band  kinetics  parameters  of  DU. 
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SURFACE  DENSITY  OF  SHEAR  BANDS  WITH  LENGTH 


IV  DETERMINATION  OF  SHEAR  BAND  KINETICS 
PARAMETERS  FROM  CFC  EXPERIMENTS 
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The  method  for  determining  the  shear  band  kinetics  parameters  from 
CFC  experiments  Is  discussed  In  this  section.  This  method,  which 
involves  fitting  the  revised  SHEAR3  model  (which  is  fully  derived  and 
described  in  Volume  II  of  this  report)  to  CFC  data,  will  be  illustrated 
here  for  both  RHA  and  4340  steel  (Rc40)  material.  First,  the  general 
procedure  is  outlined,  and  the  experimental  data  are  examined.  Then 
guidelines  for  matching  the  model  to  the  data  are  discussed,  and  the 
resulting  parameters  are  presented. 


CFC  Data 

To  determine  the  shear  band  kinetics  parameters  for  both  RHA  and 
4340  steel,  we  used  data  from  previously  performed  CFC  experiments 
spanning  a  wide  range  of  damage  from  incipient  shear  banding  to  substan¬ 
tial  fragmentation.  These  include,  for  RHA,  CFC  tests  T-2,  T-9,  and 

T-IO,  described  in  Section  III  (see  Figure  9  and  Table  2);  and  for  4)40 

14 

steel,  CFC  tests  2  and  8,  whose  geometry  and  explosive  parameters  are 
shown  in  Figure  24. 

The  primary  data  obtained  from  the  CFC  tests  are  the  final  shape  of 
the  recovered  specimen  tubes  and  the  shear  band  size  distributions. 
Preliminary  C-HEMP  simulations  are  first  performed  for  each  experiment, 
with  the  SHEAR3  model  representing  the  test  material  but  with  no  shear 
banding  allowed.  For  the  low-damage  experiments  (RHA  test  T-2  and 
steel  test  2),  the  computed  final  shape  of  the  tube  is  compared  with 


Although  shear  band  k 
for  -»340  steel  (Rc-tU) 
chat  the  changes  made 
nucleation  algorithm, 


inetics  parameters  had  been  previously  obtained*"1 
for  an  earLier  version  of  SHEAR),  it  was  relt 
in  the  revised  SHEAR),  particularly  in  the 
warranted  a  redetermination  of  these  parameters. 


51 


S' 

V 

^‘^l*  ^**  W-  ‘  -A  W.  W  _  V.  Wl  W-  **-  **-  ■*.. »*.  f  ff-  'w-'»'.~.V -  ^iitVr 


•  .  V. 
'W-V.VO.--~w 


.  -f  V. 

-V-’W-V- 


y.y-V^vjwCviVL 
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thac  obtained  experimentally  to  verify  that  the  calculation  simulates 
the  detonation  and  the  homogeneous  deformations  of  the  specimen  tube 
with  acceptable  accuracy  (such  a  comparison  is  not  useful  for  the 
higher-damage  experiments,  since  the  damage  can  significantly  affect  the 
final  shape). 

From  the  preliminary  simulations  we  obtained  histories  of  tne 
strains  in  each  damage  mode,  as  well  as  the  equivalent  plastic  strain 
and  the  various  stresses.  Plots  of  the  Mode  1  strain  histories  ' as 
defined  in  Figure  8)  are  shown  in  Figures  25  to  29  for  five  CFC  tests 
considered  here.  The  three  histories  shown  for  each  test  were  obtained 
for  ceils  in  the  midthickness  position  radially  and  at  axial  positions 
approximately  one-fourth,  one-half,  and  three-fourths  of  the  length  of 
the  tube,  respectively,  from  the  initiating  end. 

Piastic  strain  rate  histories  were  also  obtained;  a  sample  is  shown 
in  Figure  30  for  RHA  test  T-9.  The  strain  rate  shows  an  initial  pulse 
of  15-  to  20-us  duration  corresponding  to  the  passage  of  the  detonation 
wave  In  the  expLoslve.  The  second  pulse  of  significantly  Lower  peak, 
strain  rate  but  longer  pulse  duration  is  caused  by  the  rarefaction  wave 
that  is  reflected  from  the  radial  free  surface  of  the  containment 
vessel.  The  plot  in  Figure  31  of  the  strain  rate  versus  strain  also 
-.hows  this  two-pulse  structure.  The  threshold  strain  for  nucleating 
shear  bands  is  usually  reached  after  the  end  of  the  first  pulse;  there¬ 
fore,  the  lower  strain  rates  in  the  second  pulse  govern  the  nucieation 
and  growth  processes. 

Measured  shear  band  size  distributions  for  RHA  tests  T-9  and  T-1J 
have  been  shown  previously  in  Figures  lb  and  17,  respectively  (those  for 
9340  steeL  tests  2  and  8  are  reported  in  reference  14).  The  distribu¬ 
tions  for  zones  corresponding  to  the  axial  positions  for  which  the 
strain  histories  shown  above  were  obtained  were  fitted  to  the  foil  awing 
exponential  form: 


An  accurate  comparison  of  shapes  is  a  necessary  but  not  sutticient 
proof  that  the  experiment  is  well-represented  by  the  simulation. 
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FIGURE  25  CALCULATED  PLASTIC  STRAIN  HISTORIES  FOR  RHA  TEST  T-2 
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FIGURE  29  CALCULATED  PLASTIC  STRAIN  HISTORIES  FROM  4340  STEEL  TEST  8 


Ng  -  Nq  exp(-R/R1 ) , 


where  N  Is  Che  number  of  shear  bands  per  unit  volume  with  radii  greater 

O 

Chan  R,  NQ  is  Che  tncercepc  ac  zero  radius,  and  R^  is  a  shape  parameter 
for  the  disCribuCion.  NQ  and  R ^  are  decermined  in  the  fitting  process. 
These  numbers,  which  thus  characterize  the  damage  in  the  CFC  experiments, 


are  shown  in  the  "Data"  columns  in  Table  5  for  both  RHA  and  4340  steel. 


DeterminaC ion  of  Parameters 

There  are  13  parameters  provided  for  describing  shear  banding 
kinetics  in  SHEAR3.  All  must  be  specified,  but  not  all  are  critical  for 
controlling  the  damage  development.  These  parameters  are  listed  and 
briefly  defined  in  Table  6.  Some  nominal  values  for  steels  are  also 
given.  (A  more  complete  description  is  in  Volume  II  of  this  report.  ) 

In  the  following  two  subsections,  we  will  explain  how  numeric  values  are 
determined  for  these  parameters  from  CFC  data. 

The  first  two  parameters,  and  B2,  govern  growth  and  nucleation 

and  must  be  decermined  by  fitting  the  number  and  size  data  from  cylinder 
experiments,  accounting  for  the  imposed  strain  histories.  The  fitting 
procedure  is  described  in  the  next  subsection.  B-j,  the  nucleation  size 
parameter  (Rn),  is  estimated  by  extrapolating  Rp  on  a  graph  of  R[ 
versus  plastic  strain,  back  to  the  nucleation  threshold  strain  ( B ^ ,  to 
be  described).  For  the  RHA  and  4340  steel  data  considered  here  (as 
shown  in  the  "Data''  columns  in  Table  5),  it  was  difficult  to  make  an 
accurate  extrapolation  because  of  the  relative  sparseness  of  data  in  the 
intermediate  damage  range.  However,  83  =*  0.01  cm  appears  to  fit  the 
data  well. 

B^  is  the  thermal  diffusivity,  which  equals  the  thermal  conductivity 
divided  by  the  density  and  the  specific  heat  at  constant  pressure.  The 
nominal  value  for  steel  is  0.14o/s.  83  is  the  nucleation  thresnoli 
strain.  It  is  first  estimated  by  comparing  the  level  of  damage  and  peas 
plastic  strains  attained  for  a  series  of  CFC  tests  (for  example,  compar¬ 
ing  damage  in  TabLe  5  and  the  strains  in  Figures  25  through  27  tor  RHA). 
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Table  5 

SHEAR  BAND  DAMAGE  CHARACTERIZATION  PARAMETERS  FOR  RHA 
AND  4340  STEEL  FROM  CFC  DATA  AND  C-HEMP  SIMULATIONS 


Material 


Axial 

Position 


Only  incipient 
shear  band 
damage  present 


130  0.038 

290  0.047 

none 


Simulations 
N  R, 


none 

none 

0.3  0.013 


118  0.04b 

193  0.044 

none 


“5  V 

RHA 

T-10 

2.0 

4.75 

7.5 

120 

270 

5 

0. 

0. 

0. 

»*■  • 

i 

4340 

2 

5.6 

3 

0. 

pi 

Steel 

9.4 

60 

0. 

(Rc40) 

11.9 

2 

0. 

3 


Table  6 


SHEAR  BANDING  KINETICS  PARAMETERS  IN  SHEAR 3 


ParameCer  Values 


0.01  cm 


0.146  s 


0.04  cm 


Definition 


Growth  coefficient  in  equation  dR/dt  = 

B  Rdt/dt,  where  R,  £,  and  t  are  shear 
band  radius,  time,  and  plastic  strain, 
respectively 


Nucleation  coefficent 


Nucleation  size  parameter,  Rn 

Thermal  diffusivity 

Nucleation  threshold  strain 

Relative  maximum  shear  distortion 
associated  with  a  band 

Maximum  shear  band  radius  at  nucleation 


Fragment  shape  parameter  in  equation  V 
Bg(Rf)^,  where  Vv  and  are  Che  void 
volume  and  fragment  radius,  respectively 

Tangent  of  the  internal  friction  angle 
used  in  computing  shear  resistance  on 
shear  band  slip  planes 

Plastic  strain  storage  indicator 

Damage  indicator  for  triggering  slide 
lines 

Currently  unused 

Number  of  radii  used  for  each  shear  band 
size  distribution 


Nondimens ional  unless  otherwise  specified. 


The  fitting  procedure  described  below  is  then  used  to  get  a  more 
accurate  value. 

By  is  the  relative  maximum  shear  distortion  associated  with  a  band. 
It  is  determined  by  sectioning  a  recovered  CFC  tube  perpendicular  to  its 
axis  (for  Mode  L  bands)  and  measuring  the  ratio  of  the  peak  shear  dis¬ 
placement  and  the  peak  band  depth  in  the  slip  plane  (B/d  in  Figure  a). 
The  nominal  value  of  0.07  was  taken  from  previous  measurements  on  HF-1 
and  4340  steels. 

B -j  is  the  maximum  shear  band  size  at  nucleation.  This  parameter 
may  be  determined  by  observing  the  largest  bands  at  very  small  amounts 
of  strain.  The  primary  effect  of  this  parameter  is  to  determine  the 
resolution  in  the  size  distribution  for  the  calculations;  therefore,  3y 
should  be  several  times  as  large  as  Rn-  Here  we  chose  a  value  of  0.04. 

The  fragment  shape  parameter.  By,  was  chosen  to  be  1.0  to  indicate 
relatively  bulky  or  equiaxed  fragments.  Bg  is  the  tangent  of  the  angle 
of  internal  friction;  in  the  absence  of  other  data,  an  angle  of  30°  was 
used.  The  remaining  parameters,  Bjq  through  were  cnosen  to  maxe 

all  planes  active,  and  to  give  a  reasonably  fine  discretization  of  the 
size  distribution. 


Fitting  Procedure 


In  this  section,  we  discuss  the  method  we  used  to  evaluate  the  most 
critical  shear  band  kinetics  parameters  (Bp  B2,  and  By)  to  provide  the 
best  match  to  the  observed  damage.  To  perform  the  simulations  of  shear 
banding  under  an  arbitrary  Imposed  strain  histories,  we  developed  a 
small  computer  program,  NUG.EAT.  This  program  is  basically  a  simplifi¬ 
cation  of  SHEAR3.  NUCLEAT  contains  input  procedures  for  material  and 
shear  band  properties,  nucleation  and  growth  processes  for  a  single  band 
orientation,  and  the  logic  to  follow  a  prescribed  strain  path. 


To  perform  a  simulation  for  a  single  axial  position  in  one  >r 
CFC  specimen  tubes,  we  estimated  i  complete  set  of  shear  banding 
parameters  and  inserted  the  strain  history  computed  from  i  no-dam tg 
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C-HEMP  calculation  (such  as  those  in  Figures  25  through  29).  The 
initial  NUG.EAT  simulation  yielded  a  set  of  Nq  -  values  appropriate 
to  these  parameters.  Based  on  these  results  we  could  estimate  a  new  set 
of  parameters  and  repeat  the  calculation.  This  iterative  process  was 
accelerated  by  simultaneously  performing  NUCLEAT  simulations  for  all  n 
the  relevant  positions  in  all  of  the  CFC  experiments  on  a  particular 
material  (e.g. ,  all  six  4340  steel  points  or  all  nine  RHA  points 
considered  here). 

The  end  result  of  the  iterative  process  using  NUCLEAT  was  a  set  of 
shear  banding  kinetics  parameters  that  simultaneously  gave  a  best  matcn 
to  all  data  sets  for  a  single  material.  The  parameters  obtained  tor  RHA 
and  4340  steel  are  listed  in  the  row  labeled  SHEAR3  in  Table  7  (we  have 
also  included  in  this  table  the  EP  model  parameters  that  govern  the 
homogeneous  plastic  deformation  of  the  material  before  and  during  snear 
banding).  The  NQ  and  R^  values  obtained  from  the  NUCLEAT  calculations 
using  these  parameters  are  shown  under  "Simulations''  in  Table  5.  The 
correspondence  between  the  measured  and  simulated  NQ  and  Rj,  values  is 
certainly  not  precise,  but  probably  within  material  variability  (the 
variability  of  material  damage  values  is  expected  to  be  much  higher  than 
that  for  yield  strength,  moduli,  and  other  macro  parameters).  We  feel 
that  the  agreement  is  good  enough  to  indicate  that  the  growth  and 
nucleation  functions  have  forms  that  reasonably  well  represent  the 
actual  processes. 

The  next  step  in  generating  and  verifying  a  set  of  shear  banding 
kinetics  parameters  for  a  material  is  to  perform  complete  simulations  of 
the  CFC  tests,  using  C-HEMP  with  the  SHEARJ  parameters  obtained  from  the 
NUCLEAT  calculations.  The  computed  damage  from  such  simulations  would 
be  expected  to  differ  from  the  NUCLEAT  results,  especially  for  oases  n 
high  damage,  because  the  developing  damage  calculated  by  SHE.AR3  afreet 
the  stresses  and  therefore  the  subsequent  strains  in  the  C-HEMP 
simulation.  (In  contrast,  in  the  NUCLEAT  calculations,  the  impose  i 
strain  hist  ories  were  based  on  no-damage  or  Low-damage  simui it  ions. 
Unfortunately,  we  had  insufficient  time  and  funds  to  pert  arm  these 
simulations,  so  the  listed  parameters  must  be  considered  preliminarv. 
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Table  7 


SHEAR3  AND  EP  MODEL  PARAMETERS 


Parameter 

4340  Steel 

Model  Type 

R  4Q 

RHA 

SHEAR3  Shear 

Bl 

80. 0 

20.0 

banding 

b2 

1.0 

2.0 

kinetics 

b3 

0.01 

0.01 

parameters : 

S4 

0 .  14o 

0 .  I4r> 

B5 

e 

• 

ro 

•F 

0.3  5 

B6 

0.07 

0.07 

B7 

0.04 

0.04 

B8 

1.0 

1.0 

b9 

0.577 

0.577 

B10 

2.0 

2.0 

8  L 1 

2.0 

2.0 

S12 

0.0 

0.0 

B1 3 

8.0 

8.0 

EP  Yield 

10.  3 

b .  8 

Strengch 

10.7 

7.  7 

(Kbar ) 

13.7 

8.  7 

9.5 

10.  -4 

11.  1 

12.  -4 

14.0 

V  COMPUTER  SIMULATION  OF  ARMOR  PENETRATION 


A  series  of  armor  penetration  simulations,  involving  impact 
velocities  both  near  and  well  above  the  ballistic  limit,  were  performed 
using  the  newly  developed  C-HEMP  code.  There  were  three  primary 
purposes  for  these  calculations.  The  first  was  to  test  some  of  the 
capabilities  of  C-HEMP  (in  particular  the  slide  line  and  rezoning 
features)  under  the  extreme  deformation  conditions  that  exist  in  a 
penetration  scenario.  The  second  was  Co  study  the  effect  of  the 
inclusion  of  the  SHEAR3  damage  model  upon  the  penetration  simulations. 
And  the  third  was  to  compare,  at  least  qualitatively,  the  simulated 
results  with  the  results  of  previously  performed  experiments  having  a 
similar  geometry  and  impact  velocity. 

All  of  the  simulations  involved  Che  normal  impact  of  a  cylindrical 
4 340  steel  (Rc40)  penetrator  into  a  circular  RHA  plate,  a  scenario  for 
which  we  have  a  large  experimental  data  base.  The  geometry  and  impact 
velocity  (V^)  for  the  five  simulations  that  will  be  discussed  in  this 
section  are  shown  in  Table  8  and  Figure  32.  Simulations  No.  1  and  2 
represent  che  case  of  a  long-rod  penetrator  impacting  at  well  above  the 

ballistic  limit.  The  corresponding  experiment  was  performed  at  the  HRL 

1 

quarter-scale  facility  during  the  second  year  of  this  program  and  is 
depicted  in  Figure  33.  Simulations  No.  3,  4,  and  5  represent  the  case 
of  a  shorter  penetrator  near  or  below  the  ballistic  limit.  The  corres¬ 
ponding  experiments  were  performed  ac  SRI  International  during  the  first 
year  of  this  program,^  and  the  results  are  depicted  in  Figure  34. 

The  cylindrical  symmetry  of  the  normal  penetration  scenario  allows 
the  two-dimensional  C-HEMP  code  co  be  used  in  its  axisymmetrie  mode.  \ 
sample  initial  cell  layout  is  shown  in  Figure  35.  Each  quadri later  i  1 
ceLl  actually  represents  the  toroidal  volume  obtained  by  rotating  thit 
:e  1 1  iriund  Che  centerline  'axis  of  symmetry!.  A  s  I  i  it*  line  extends 
i  1  >ng  the  impact  interface,  allowing  the  end  jf  the  penetrat  »r  to  s ;  .  :*.• 
wirnout  fricti>n  along  the  front  tace  of  the  irmor, 

0/ 
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FIGURE  35 
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PORTION  OF  ORIGINAL  CELL  LAYOUT 
FOR  PENETRATION  SIMULATIONS 


The  cells  in  the  vicinity  of  the  impact  interface  are  subject  to 
severe  shear  deformations,  and  because  the  time  step  (the  time  interval 
between  successive  computational  cycles)  must  be  less  than  the  time  for 
an  incremental  stress  wave  to  propagate  across  any  cell,  the  time  step 
decreases  precipitously,  thus  ending  the  calculation. 

We  have  used  a  three-part  strategy  to  mitigate  this  effect.  First, 
all  of  the  cells  in  the  vicinity  of  the  impact  interface  are  continuously 
rezoned.  Second,  the  cells  are  allowed  to  collapse  into  triangles;  if 
one  node  of  the  quadrilateral  is  close  to  penetrating  one  of  the 
opposite  sides,  that  node  is  held  fixed  in  the  direction  perpendicular 
to  the  side  it  is  about  to  penetrate,  and  that  node  is  eliminated  from 
the  time  step  calculation.  Third,  as  a  last  resort,  any  cell  whose 
calculated  time  step  falls  below  a  user-specified  value  is  eliminated 
from  further  calculational  cycles.  As  a  result,  subsequent  motions  and 
deformations  in  the  vicinity  of  such  a  cell  will  be  erroneous,  but  there 
should  be  little  deleterious  effect  on  regions  far  removed. 

The  SHEAR3  model  (which  is  described  in  detail  in  Volume  II  of  this 
report)  was  used  for  both  the  armor  and  penetrator  materials  in  three  of 
the  simulations  (Nos.  1,  3,  and  5)  to  calculate  the  nucleation,  growth, 
and  coalescense  of  shear  bands,  and  the  resultant  loss  of  shear  strength. 
The  shear  banding  kinetics  parameters  used  in  SHEAR3  were  obtained  from 
previous  CFC  experiments  and  computational  modeling  (as  described  in 
Section  IV),  and  are  listed  in  Table  7. 

Figure  36  shows  the  shear  band  damage  profile  for  simulation  No. 1 
(V^  *  1550  m/s  with  SHEAR3)  at  11.1  ps  after  impact,  which  is  as  far  as 
the  calculation  proceeded  before  the  time  step  became  excessively  small. 
The  projectile  has  penetrated  about  40%  of  the  way  into  the  target, 
while  the  back  surface  of  the  target  has  experienced  only  very  small 
motions.  The  kinetic  energy  of  the  projectile  is  still  approximately 
38%  of  its  original  value.  The  pointed  projectile  tip  is  probably  a 
result  of  insufficient  computational  cell  resolution  in  the  radial 
direction  in  the  penetrator.  A  simulation  with  better  resolution 
exhibited  less  of  a  pointed  tip,  but  resulted  in  earlier  time  step 
degradation. 


The  dark  shaded  regions  in  Figure  36  (and  in  subsequent  figures) 
represent  cells  that  have  been  completely  damaged  by  shear  banding  and 
thus  no  longer  retain  any  shear  (or  tensile)  strength.  The  lighter 
shaded  regions  represent  cells  that  have  undergone  partial  shear  banding 
damage.  The  damaged  region  occupies  the  impact  end  of  the  projectile  as 
well  as  a  broad,  relatively  uniform  region  of  the  target  adjacent  to  the 
impact  zone.  From  the  progression  of  damage  at  times  less  than  11  us 
and  from  the  relatively  small  loss  of  projectile  kinetic  energy  at 
11  ps,  one  can  roughly  extrapolate  the  progression  of  damage  at  later 
times.  It  appears  that  the  entire  thickness  of  the  target  plate  in  the 
vicinity  of  impact  will  become  fully  damaged,  while  the  projectile  still 
recains  a  large  fraction  of  its  kinetic  energy,  thus  reasonably  well 
simulating  (at  least  qualitatively)  the  large  volume  of  small  fragments 
observed  in  Figure  33. 

Figure  37  shows  the  damage  profiles  for  simulation  No.  3  ( V ^  =  7  50 
m/s)  12.3  ps  after  impact.  The  projectile  has  penetrated  60%  of  the  way 
into  the  target,  the  back  surface  of  the  target  has  moved  about  40%  of 
its  thickness,  and  Che  kinetic  energy  of  the  projectile  has  dropped  to 
38%  of  its  original  value. 

The  damaged  region  in  Che  projectile  is  similar  to  that  for  the 
higher  velocity  simulation.  But  in  the  target  the  situation  is  quite 
different.  Rather  than  occupying  a  broad  region  adjacent  to  the  impact 
zone,  the  fully  damaged  area  localizes  along  a  narrow  band  (one  cell 
thick)  extending  toward  the  back  of  the  plate.  The  band's  distance  from 
the  impact  axis  is  slightly  greater  than  the  projectile  radius.  From 
the  progression  of  damage  at  earlier  times  and  from  the  remaining 
projectile  kinetic  energy,  it  is  clear  that  a  continuation  of  the 
simulation  (if  the  time  step  would  have  permitted  it)  would  show  the 
damage  band  reaching  the  back  surface  of  the  target  and  forming  a  plug. 
The  plug  and  the  remainder  of  the  target  would  be  largely  undamaged  by 
shear  banding,  a  result  similar  to  that  seen  experimentally  for  a  792 
m/s  impact  [Figure  34(c)]. 


MA-7893-175 


FIGURE  37  PROFILE  OF  SIMULATION  NO.  3  (V,  =  750  m  s)  12.3  ps  AFTER  IMPACT 

Dark  shaded  regions  indicate  computational  cells  that  are  completely  damaged 
by  shear  banding.  Light  shaded  regions  indicate  cells  that  have  experienced 
some  shear  banding  failure. 
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Figure  38  shows  Che  damage  profile  for  simulation  No.  5  (V^  =  out) 
m/s)  ac  12.5  ,is  after  impact.  The  front  of  the  projectile  Is  fully 
damaged  by  shear  banding,  but  the  target  has  undergone  insufficient 
strains  to  nucleate  any  shear  bands.  The  projectile  has  penetrated  3QS 
of  the  way  into  the  target,  the  back  surface  of  the  target  has  moved 
about  1/6  of  its  thickness,  and  the  projectile  kinetic  energy  is  only 
21%  of  its  original  value.  It  is  therefore  unlikely  that  a  continuation 
of  the  calculation  would  produce  any  significant  shear  banding  in  the 
target;  the  final  profile  would  appear  similar  to  that  of  the  experi¬ 
mental  655  m/s  impact  (Figure  34(a)], 

These  lower  impact  velocity  simulations  qualitatively  predict  the 
progression  of  damage  we  have  observed  in  a  target  in  the  vicinity  of 
the  ballistic  limit — from  homogenous  plastic  deformation,  to  shear 
strain  localization,  to  complete  plugging.  In  addition,  the  calcula¬ 
tions  support  our  experimental  finding  that  the  ballistic  limit  for  this 
particular  geometry  is  between  710  and  792  m/s.*  However,  considering 
the  large  cell  size  and  the  preliminary  nature  of  the  shear  banding 
parameters  used  In  the  calculations,  it  was  not  expected  that  we  would 
accurately  predict  a  ballistic  limit.  What  is  significant  here,  we 
believe,  is  the  ability  of  the  simulations  to  correctly  predict,  using 
the  same  damage  model  and  kinetics  parameters,  the  qualitative  differ¬ 
ence  between  the  penetration  phenomenology  observed  at  the  two  widely 
separated  impact  velocity  regimes. 

To  more  fully  assess  the  impact  of  the  SHEAR3  damage  model  on  the 
simulated  penetration  at  velocities  both  near  and  well  above  the 
ballistic  limit,  we  repeated  two  simulations  replacing  SHEAR3  with  the 
elastic-plastic  (EP)  model  (Nos.  2  and  4).  The  EP  model  allows  homo¬ 
geneous  shear  deformation  along  a  work-hardening  yield  curve,  but  no 
shear  banding.  The  EP  model  parameters,  determined  from  quasi-static 


*The  ballistic  limit  is  a  statistical  quantity  and  can  be  determined 
only  from  many  tests. 
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FIGURE  38  PROFILE  OF  SIMULATION  NO.  5  (V,  =  500  m/s)  12.5  ys  AFTER  IMPACT 

Dark  shaded  regions  indicate  computational  cells  that  are  completely  damaged 
by  shear  banding.  Light  shaded  regions  indicate  cells  that  have  experienced 


some  shear  banding  failure. 


tensile  tests,  are  also  listed  in  Table  7  in  terras  of  yield  strength  as 
a  function  of  equivalent  plastic  strain.  The  constitutive  paths 
described  by  the  SHEAR.3  and  EP  models  are  identical  before  the  onset  of 
shear  banding.  Comparisons  between  simulations  using  SHEAR3  and  other¬ 
wise  identical  simulations  using  the  EP  constitutive  model  are  described 
below. 

For  the  higher  velocity  simulations  (V^  =*  1530  m/s),  there  appears 
to  be  no  significant  difference  in  the  cell  profiles  at  about  11.5  -s 
after  impact  (see  Figure  39).  Simulation  No.  1,  using  SHEAR3,  resulted 
in  slightly  greater  deformation  of  both  the  projectile  and  target  in  the 
region  near  the  impact  interface,  and  slightly  less  penetration,  than 
simulation  No. 2,  using  EP. 

There  was  no  significant  difference  in  the  simulated  deformations 
occurring  early  in  a  high-velocity  penetration  scenario  between  models 
having  different  yield  and  shear  failure  algorithms.  This  is  not 
surprising  because  Che  early  stresses  are  very  high  compared  with  the 
yield  strengths,  and  thus  the  compressibility  and  momentum  relations 
largely  govern  the  deformations.  Of  course,  only  SHEAR3  can  model  the 
actual  shear  band  failure  processes  and  the  resultant  reduction  in  shear 
strength,  and  therefore  be  able  to  predict  deformations  later  in  the 
penetration  (when  material  strengths  play  a  more  important  role),  as 
well  as  fragmentation. 

For  the  lower  velocity  simulations  (Vi  =  750  m/s),  large 
differences  in  deformation  do  occur  in  the  first  12  ps  after  impact  (see 
Figure  40).  Simulation  No.  3,  using  SHEAR3,  resulted  in  far  greater 
deformations  in  the  impact  region  for  both  the  projectile  and  target 
than  simulation  No.  4,  which  used  EP.  More  significantly,  the  simula¬ 
tion  using  EP  did  not  experience  any  of  the  localization  of  shear  damage 
discussed  above  (and  shown  in  Figure  37),  which  would  eventually  lead  to 
a  plugging  type  of  target  failure.  As  a  result,  the  projectile  penetra¬ 
tion  and  the  target  back  surface  deflection  are  about  25%  smaller  in 
simulation  No.  4  than  in  simulation  No.  3. 


Upper  Profile  —  Simulation  No.  4  (with  EP  model),  12.5  us  after  Impact 


Lower  Profile  —  Simulation  No.  3  (with  SHEAR3  model),  12.3  us  after  Impact 

MA-7893-197 

FIGURE  40  EFFECT  OF  SHEAR  BANDING  DAMAGE  ON  SIMULATED  PENETRATION 
AT  AN  IMPACT  VELOCITY  OF  750  m/s 


In  conclusion,  although  we  have  not  been  able  to  simulate  to 
completion  an  armor  penetration  by  projectile  impact,  we  have,  by  use  of 
the  C-HEMP  code  with  the  SHEAR3  model,  extended  the  simulation  suf¬ 
ficiently  to  demonstrate  that  it  can  predict,  at  least  qualitatively, 
much  of  the  armor  penetration  phenomenology  that  we  have  observed 
experimentally. 
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VI  RECOMMENDATIONS  FOR  FUTURE  JORX 


As  documented  in  this  report,  SRI  International  has  made  u^niri  :jr. 
progress  in  developing  a  phenomenologically  sound,  mater ial-pr jpert- -Das 
computational  method  for  predicting  armor  penetration  by  projectile 
impact.  At  this  point,  the  method  can  be  used  to  simulate  moueritel.' 
deep  penetrations  caused  by  normal  impact  of  long-rod  projectiles,  mu 
to  correctly  predict  the  mode  of  penetration  as  well  as  rough  estimates 
of  the  resultant  fragment  size  and  velocity  distributions.  However, 
much  needs  to  be  done  to  extend  the  method  to  deep  penetrations  and 
oblique  impacts,  and  to  make  it  readily  useable  by  DoD  personnel  to 
facilitate  the  design  of  new  armor,  penetrator,  or  fragmentation 
systems.  In  this  chapter  we  present  our  recommendations  for  achieving 
this  goal. 

The  basic  aim  of  the  future  work  we  are  recommending  is  to  improve, 
simplify,  and  transfer  to  BRL  an  improved  capability  for  computationally 
predicting  armor  penetration  and  the  resultant  fragment  size,  shape,  and 
velocity  distributions.  The  approach  we  suggest  is  as  follows: 

•  Determine  what  modifications  can  be  made  to  the  shear 
band  model,  the  computer  code  (C-HEMP),  and  the  experi¬ 
mental  procedures  for  characterizing  a  material,  to 
Improve  accuracy,  facilitate  use,  and  extend  the 
capability  to  deep  penetrations  and  oblique  impacts. 

•  Implement  these  modifications,  testing  their  effective¬ 
ness  by  computational  simulations  of  previously  performed 
penetration  experiments. 

•  Apply  the  computational  method,  with  the  modifications, 
to  the  predictions  of  various  penetration  (or  ocher 
appropriate)  scenarios  of  interest  to  BRL. 

•  Document  the  model,  the  code,  and  the  procedures  for 
material  characterization,  and  transfer  this  capability 
to  BRL  In  a  manner  that  will  make  it  readily  available 
for  use  In  design  calculations. 


*>  *. 


All  modif ications  to  Che  computational  method  would  remain  within 
our  overaLl  philosophy  of  basing  che  model  on  a  quantitative  description 
of  the  phenomenological  processes  Chat  do  occur  in  an  actual  penetration. 
We  can  thus  be  assured  chat  Che  model  will  be  applicable  to  arbitrary 
geometries. 

We  have  divided  the  recommended  future  work,  into  six  suggested 
tasks,  which  are  described  below. 

Task  1:  Refinements  to  Shear  Band  Model 

Examination  of  experiments  Involving  both  normal  and  oblique  impact 
of  long-rod  penetrators  into  RHA  plates^  have  revealed  several  phenom¬ 
enological  features  Chat  are  either  not  yet  incorporated  in  our  shear 
banding  model,  or  if  incorporated,  have  not  been  fully  tested  or  exer¬ 
cised.  These  features  include  the  significance  of  shear  band  anisotropy 
to  che  penetration  and  fragmentation  process,  tensile  separation  along 
an  already-formed  shear  band,  and  the  possibility  that  a  specific  point 
in  a  material  may  begin  to  fail  by  tension  before  it  has  experienced 
sufficient  shearing  to  form  shear  bands.  The  last  feature  is  extremely 
important  in  fragmenting  rounds,  but  may  also  be  important  in  penetra¬ 
tion  scenarios  In  which  tensile  failure  in  the  target  precedes  plug 
formation. 

SRI  International  has  recently  been  developing,  primarily  under 

★ 

another  program,  a  model  for  shear  banding  in  an  anisotropic  material, 
SHEARS,  which  describes  the  nucleation  and  growth  of  bands  in  as  many  as 
seven  distinct  orientation  modes.  This  model  is  anisotropic:  It 
relaxes  shear  stresses  and  reduces  shear  strengths  only  on  the  specific 
planes  that  are  undergoing  shear  banding,  and  it  allows  for  different 
nucleation  and  growth  parameters  for  different  planes  (to  account  r  r 
material  anisotropy).  It  also  aJ lows  the  shear  bands  of  a  pirt. 
orientation  mode  to  change  direction  as  their  cell  der  >r  < . 

*Array  Research  Office  Contract  OAAO-d--'. 


of  Che  bands  in  Che  ocher  orlenCaclon  modes.  However,  chese  anisotropic 
capabilides  have  noc  been  fully  cesCed  or  exercised. 

There  exiscs,  from  che  current  program  (see  Section  III),  CFC  data 
on  RHA  tubes  whose  axes  have  been  cut  parallel  to  and  perpendicular  to 
the  rolling  plane.  Up  to  now,  we  have  only  considered  orientation  Mode 
1  bands  when  attempting  to  computationally  simulate  the  experiments.  We 
now  recommend  testing  the  ability  of  the  code  to  predict  shear  banding 
in  other  than  Mode  1  orientations,  using  the  pararaenters  determined  from 
Mode  1  simulations  in  the  two  directions.  This  would  provide  a  fuller 
check  on  the  anisotropic  capabilities  of  SHEAR4. 

Furthermore,  in  all  previous  simulations  using  either  of  our  two 
shear  banding  models,  a  specific  computational  cell  is  allowed  to 
nucleate  and  grow  shear  bands  but  not  tensile  voids,  even  if  the  stress 
and  strain  history  for  that  cell  warrants  initial  tensile  failure.  For 
scenarios  that  include  simultaneously  both  shear  banding  and  tensile 
failure  (e.g.,  fragmenting  rounds),  the  programmer  must  specify  in 
advance  which  cells  will  follow  the  shear  banding  model  and  which  will 
follow  a  tensile  failure  model.  We  therefore  recommend  examining  the 
possibility  of  having  each  cell  determine,  from  its  loading  history  and 
material  failure  properties,  what  type  of  failure  it  undergoes,  and 
switch  from  one  type  to  another  if  the  loading  changes  appropriately. 

If  this  feature  appears  to  be  feasible,  we  recommend  implementing  it  in 
the  SHEAR4  model  and  testing  it  against  existing  experimental  data. 

From  phenomenological  experiments,  we  have  seen  that  armor  and 
penetrator  fragmentation  is,  in  most  cases,  caused  by  the  intersection 
of  shear  bands  propagating  in  different  directions.  In  our  current 
shear  banding  models,  fragmentation  in  a  particular  computational  cell 
is  said  to  have  occurred  when  a  damage  parameter  (a  simple  function  of 
the  density  and  size  of  the  shear  bands  in  that  cell)  reaches  unity. 

The  resultant  fragment  size  distribution  comes  from  the  final  shear  band 
size  distribution  before  fragmentation  (many  small  bands  produce  many 
small  fragments,  while  a  fewer  and  larger  bands  produce  fewer  and  larger 
fragments). 


However,  no  attempt  is  made  to  describe  the  shape  or  aspect  ratio 
of  these  fragments,  although  pertinent  information  is  available,  in 
terms  of  the  size  distributions  of  the  shear  bands  in  the  various 
orientation  modes  present  within  that  cell.  For  example,  if  there  are 
similar  distributions  of  bands  in  several  different  orientations,  then 
chunky  fragments  (aspect  ratio  nearly  unity)  might  be  expected.  But  if 
one  orientation  had  a  much  lower  density  of  bands,  then  elongated  frag¬ 
ments  (large  aspect  ratios)  might  be  expected.  We  therefore  recommend 
exploring  methods  of  improving  the  simple  damage  parameter  determination 
of  cell  fragmentation  and  implementing  a  fragment  shape  prediction 
capability  in  the  SHEAR4  model. 


Task  2:  Improvements  to  C-HEMP  for  Armor  Penetration  Calculations 


The  current  version  of  C-HEMP  incorporates  many  of  the  features 
considered  important  for  armor  penetration  calculations.  The  finite 
element  cell  numbering  scheme,  which  makes  cell  numbering  independent  of 
cell  location,  allows  more  flexibility  for  slide  line  insertion  or 
separating  (and  eliminating)  spalled  regions.  The  capability  to  handle 
complex  material  models  makes  it  possible  to  use  SHEAR3,  SHEAR4,  and 
other  MSFM  models,  such  as  the  NAG  ductile  and  brittle  failure  models, 
as  well  as  any  new  models  to  be  developed  in  the  future.  The  current 
slide  line  capability  allows  slip  (with  or  without  a  frictional 
component)  along  user-located  slide  lines  that  are  either  already 
"unzipped"  or  "unzip”  as  the  computation  proceeds.  And  the  automatic 
rezoner  compensates,  to  a  certain  extent,  for  the  large  cell  distortions 
involved  in  a  penetration  simulation. 


However,  the  main  problem  with  C-HEMP  is  its  inability  to  complete 
deep  penetrations  simulations  because  extreme  cell  distortions  decrease 
the  allowable  time  step  to  unworkably  small  values.  The  usual  solution 
to  this  problem  is  to  replace  quadrilateral  cells  with  triangular  celLs. 
While  the  triangular  cells'  resistance  to  collapse  (as  compared  to  a 
quadrilateral  cell,  in  which  one  node  may  collapse  onto  an  opposite 
side)  may  in  some  situations  cause  unrealistically  high  pressures,  it 


keeps  Che  cime  seep  reasonably  large  and  thus  allows  Che  calculations  co 
proceed  Co  conclusion.  We  expect  that  the  excessive  pressures  will  be 
less  of  a  factor  in  Che  primarily  shear-strain-driven  shear  banding 
models  chan  in  the  primarily  pressure-driven  pore  compaction  model,  for 
example. 

What  appears  to  be  most  promising  is  a  mixed  approach  whereby  both 
triangular  and  quadrilateral  cells  are  used  in  the  same  calculation: 
the  triangular  cells  are  used  in  the  regions  of  expected  extreme  defor¬ 
mations.  Such  an  approach  was  used  for  a  deep  penetration  simulation 

★ 

with  the  SWE2D  code,  as  shown  in  Figures  41  and  42.  Of  course,  this 
code  does  noC  Include  SHEAR3  and  thus  cannot  predict  plugging  or  frag¬ 
mentation,  but  it  can  handle  extreme  cell  deformations  (similar  results 
have  been  obtained  with  other  triangular  cell  codes,  such  as  BRL'S  EPIC 
code).  Combining  a  triangular  cell  capability  with  the  shear  banding 
models  and  other  features  currently  in  C-HEMP  would  be  the  next  logical 
step  toward  completing  a  computational  method  for  predicting  deep 
penetrations. 

In  addition  to  Che  mixed  triangular  and  quadrilateral  cell  feature, 
a  more  complex  slide  line  capability  is  needed  to  better  simulate  both 
normal  and  oblique  impacts.  For  normal  impacts  of  long-rod  penetrators 
at  velocities  somewhat  greater  than  the  ballistic  limit  of  the  armor 
plate,  penetration  occurs  by  the  propagation  of  a  shear  band  across  the 
thickness  of  the  armor  plate  to  form  a  plug.  When  this  band  reaches  the 
back  surface  of  the  armor,  the  plug  is  removed. 

One  method  of  computationally  simulating  this  scenario  is  to  insert 
a  slide  line  at  an  appropriate  position  between  two  rows  of  armor  plate 
cells  parallel  to  the  direction  of  impact.  This  slide  line  is  initially 
inactive,  but  is  gradually  activated  along  its  length  by  a  criterion 
related  co  chat  of  nucleation  and  growth  in  our  shear  band  model. 


Written  by  Thomas  Cooper  at  the  Swedish  Detonic  Institute. 
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FIGURE  41  INITIAL  CELL  PROFILES  FOR  SWE2D  DEEP  PENETRATION 
SIMULATION.  USING  COMBINATION  OF  TRIANGULAR 
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AND  DETAIL  (b)  OF  HIGHLY  DEFORMED  REGION 


Another  level  of  slide  line  complexity  arises  because  the  location 
of  the  unzipping  slide  line  cannot  necessarily  be  determined  in  advance 
of  the  computation.  If  a  blunt-nosed  hardened  penetrator  impacts 
normally  into  a  4340  steel  (Rc40)  armor  plate,  for  example,  a  shear  band 
nucleates  in  the  armor  at  the  very  edge  of  the  penetrator,  forming  a 
plug  with  the  same  diameter  as  that  of  the  penetrator.  For  this  case, 
the  slide  line  can  be  user-located  in  advance.  However,  for  a  round- 
nosed  penetrator  or  for  a  more  realistic  armor  material,  such  as  RHA, 
significant  material  distortion  may  occur  in  both  the  armor  and 
penetrator  before  the  nucleation  of  a  shear  band  in  the  armor,  and  plugs 
have  been  formed  with  a  diameter  significantly  greater  than  that  of  the 
penetrator. ^ 

It  would  thus  be  desirable  to  have  the  code  select  the  position  of 
the  slide  line  as  the  calculation  progresses.  This  "self-locating" 
slide  line  could  be  positioned  by  the  code  adjacent  to  the  first  ceil 
that  satisfies  a  shear  banding  damage  criterion  that  corresponds  to 
nearly  complete  loss  of  shear  strength  in  the  direction  of  the  slide 
line.  Implementing  this  feature  should  be  relatively  straightforward, 
since  the  finite  element  numbering  scheme  would  allow  for  creating  extra 
cell  nodes  as  required  by  adding  a  slide  line  in  the  middle  of  a 
calculation.  Once  the  slide  line  is  positioned,  it  would  proceed  to 
unzip  from  node  to  node  as  the  cells  adjacent  to  the  nodes  satisfy  the 
same  damage  criterion. 

Not  only  is  the  self-locating  unzipping  slide  line  more  realistic 
from  a  phenomenological  viewpoint  and  requires  less  user  intervention, 
but  it  has  a  third  advantage.  As  described  in  Section  V,  simulations  of 
normal  rod  impacts  without  the  use  of  slide  lines  have  resulted  in  the 
localization  of  shear  banding  damage  in  a  particular  row  of  cells 
parallel  to  the  impact  direction.  The  ensuing  loss  of  shear  strength 
causes  these  cells  to  undergo  extreme  shear  distortions,  significant Lv 
reducing  the  time  step  between  cycles  and  thereby  greatly  increasing  the 


calculation  time.  A  self-* locating  slide  line  would  alleviate  this 
problem,  and  allow  for  complete  simulation  by  letting  the  plug  separate 
from  the  target. 


■  We  have  thus  far  assumed  that  the  lateral  surfaces  of  the  plug  are 

|  relatively  straight  and  perpendicular  to  the  armor  plate.  Such  is  not 

|  '  the  case  in  many  armor  penetration  scenarios.  During  normal  impact  at 

I  velocities  below  or  slightly  above  the  ballistic  limit,  the  shear  bands 

i 

1  that  form  the  lateral  edges  of  the  plugs  turn  inward  (toward  the 

'  29 

|  penetracor  axis)  as  they  grow  toward  the  back  of  the  armor.  And 

during  oblique  impact,  the  penetrator  may,  as  it  proceeds  through  the 
'  target  plate,  rotate  to  an  orientation  more  normal  to  the  target  face. 

So  the  hole  formed  by  the  plug,  or  the  fragmented  material  removed,  has 
surfaces  with  curvatures  that  are  certainly  not  predictable. 

We  are  thus  led  to  a  further  level  of  slide  line  complexity — a 
slide  line  that  changes  its  direction  as  it  propagates,  according  to  the 
stress  and  strain  fields  present.  This  “self-directed"  (or  "wandering”) 
slide  line  feature  will  not  be  simple  to  implement.  One  problem  is 
handling  the  changing  cell  layout  as  the  slide  line  propagates  through 
the  grid,  cutting  through  the  middle  of  cells  instead  of  merely  separat¬ 
ing  adjacent  cells.  In  principle,  the  finite  element  numbering  scheme 
can  handle  such  a  procedure  by  creating  new  nodes  wherever  they  are 
required  (refer  to  Appendix  F  of  reference  16),  but  significant  new 
programing  will  be  required.  Another  problem,  and  perhaps  the  most 
difficult,  will  be  to  determine  the  criterion  that  controls  the  direc¬ 
tion  the  slide  line  takes.  The  direction  of  the  plane  of  maximum  shear 
strain  would  be  one  ingredient  in  this  criterion,  but  it  is  clear  that 
the  physics  of  the  problem  will  have  to  be  studied  more. 

So  far,  we  have  discussed  only  nonintersecting  slide  lines  (these 
are  often  referred  to  as  "single"  slide  lines,  but  there  may  be  more 
than  one  in  a  calculation,  provided  they  don't  intersect).  But  for  a 
penetration  scenario  that  includes  a  combination  of  penetrator  deforma¬ 
tion  and  removal  of  a  plug  (or  any  fragmented  region) — and  this  includes 
almost  all  realistic  penetrations — there  is  a  need  for  intersecting  or 
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'double'*  slide  lines.  One  slide  line,  located  within  the  armor  plate, 
allows  the  plug  to  slip  out  from  the  rest  of  the  plate.  The  other  slide 
line,  whose  slip  surfaces  initially  comprise  the  front  surface  of  the 
armor  plate  and  the  impact  and  lateral  edges  of  the  penetrator,  respec¬ 
tively,  allows  the  distorting  penetrator  to  slide  with  respect  to  the 
distorting  target.  Because  both  processes  may  occur  simultaneously,  and 
because  the  two  slide  lines  intersect,  one  pair  of  nodes  of  one  slide 
line  will  need  to  slip  along  the  other  line.  And  therein  lies  the 
complication. 

We  have  very  recently  attempted  a  double  slide  line  calculation 
with  an  unzipping  slide  line  as  one  of  two  intersecting  lines.  It 
appeared  to  work  well,  but  this  feature  needs  to  be  tested  and  exercised 
in  different  goemetries  to  complete  the  Implementation. 

Task.  3;  Proof  Tests  of  the  Computational  Method 

To  determine  the  accuracy  of  our  computational  method  for  predicting 
penetration  and  fragmentation,  we  recommend  that  proof  tests  be  performed 
by  simulating  several  penetration  scenarios  for  which  adequate  data 
exists.  These  scenarios  should  cover  a  reasonably  wide  range  of  impact 
velocity  and  resultant  damage  to  exercise  as  many  different  features  of 
the  computational  method  as  possible.  But  the  geometries  should  be  kept 
as  simple  as  possible  so  that  these  proof  tests  may  be  performed  before 
some  of  the  modifications  described  in  the  first  two  tasks  are 
implemented.  By  simulating  only  normal  impacts,  for  example,  we  can 
perform  the  proof  tests  without  first  having  to  implement  the  self- 
directing  slide  line  feature. 

Two  experiments  performed  during  the  current  BRL  program  are 
particularly  suitable  for  testing  a  wide  range  of  computational 
features.  Both  involve  the  normal  impact  of  a  olunt-nosed  4340  steeL 
rod  into  an  RHA  plate.  Thus,  4340  steel  (Rc40)  and  RHA  (cut  from  LO-cro- 
thick  plates)  are  currently  the  materials  whose  shear  banding  kinetics 
have  been  best  characterized. 


* 


One  Impact,  at  0.792  km/s,  just  above  the  ballistic  limit,  resulted 
in  the  removal  of  a  plug,  but  there  was  no  other  damage  (except  homo¬ 
geneous  plastic  deformation)  in  the  plug  or  the  remainder  of  the  armor 
plate  (see  Figure  34).  The  penetrator  remained  intact,  embedded  in  the 
plug,  with  some  shear  banding  damage  in  the  impact  end.  The  other 
impact,  at  1.55  km/s,  resulted  in  a  large  amount  of  downrange  fragmenta¬ 
tion  from  both  the  armor  and  the  penetrator,  mostly  in  the  form  of  small 
fragments  (see  Figure  33).  The  region  in  the  target  plate  adjacent  to 
the  penetration  hole  was  heavily  shear  banded. 


We  recommend  performing  proof  tests  of  C-HEMP  and  SHEAR4  by 
simulating  the  experiments  described  above  (or  other  penetration 
experiments).  It  may  be  necessary  to  complete  implementation  of  some  of 
the  modifications  described  in  the  first  two  tasks  (the  double  slide 
line  feature,  for  example)  before  performing  these  proof  tests. 
Furthermore,  as  some  of  the  more  advanced  features  become  implemented  in 
C-HEMP  and  SHEAR4,  additional  proof  tests  will  be  needed  (an  oblique 
impact,  for  example,  to  test  the  self-directing  slide  line). 


Finally,  we  recommend  one  other  type  of  proof  test.  It  relates  to 


the  use  of  a  two-dimensional  code  to  simulate  a  three-dimensional 


problem— the  oblique  impact  of  a  rod  on  a  plate.  We  plan  to  solve  this 
problem  with  the  standard  planar  solution:  the  oblique  impact  of  a 
plate  of  infinite  width  upon  another  plate  of  infinite  width.  But  we 
need  to  know  how  precise  this  approximation  is,  and  in  particular,  how 
accurate  are  the  final  calculated  results  (e.g.,  fragment  size  and 
velocity  distributions). 


Such  a  determination  can  be  made  by  comparing  the  results  of  two 
simulations  of  a  set  of  zero-obliquity  rod-tnto-plate  impacts  (such  as 
the  pair  described  above).  One  simulation  would  use  the  true  axisym- 


metric  geometry,  while  the  other  would  use  the  planar  approximation.  A 
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comparative  study  by  Zukas  using  ogival  projectiles  showed  that  the 
planar  approximation  was  quite  accurate  at  early  times,  but  that  its 
accurary  degraded  with  increasing  time  because  it  neglected  such 
Important  physical  phenomena  as  the  out-of-plane  motions  leading  to 
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laCeral  stress  relaxation.  We  recommend  performing  a  similar  compara¬ 
tive  study,  using  Che  blunt-nosed  projectile  impacts  described  above,  to 
assess  the  accuracy  of  the  planar  approximation  in  oblique  Impact 
simulations  with  SHEAR4  and  C-HEMP. 

Task  4:  Simplification  of  Computational  Procedure 

The  ultimate  goal  of  this  effort  is  to  develop  a  procedure  avail¬ 
able  for  routine  use  by  both  SRI  and  BRL  (and  eventually  other  DoD 
agencies)  for  armor  and  penetrator  design  calculations,  as  well  as  for 
other  fragmentation  applications.  We  recommend  simplifying  both  the 
calculational  models  and  the  experimental  procedures  necessary  to  char¬ 
acterize  a  material  to  be  modeled,  without  sacrificing  the  procedures 
integrity  or  accuracy.  Below  we  describe  various  methods  of  effecting 
this  simplification. 

One  method  is  to  simplify  the  computational  model  for  shear  banding 
by  eliminating  from  the  nucleation,  growth,  or  coalescence  criteria 
those  variables  chat  do  not  significantly  affect  the  outcome.  These  may 
Include  material  parameters  obtained  from  the  characterizing  experi¬ 
ments,  as  well  as  constants  obtained  from  the  literature.  Determining 
which  variables  can  be  safely  eliminated  can  be  done  by  systematically 
varying  the  input  parameters  used  in  a  series  of  penetration  simulations 
and  observing  the  effect  on  the  final  results  (such  as  the  fragment  size 
distribution).  If  changes  in  a  particular  parameter  have  little  or  no 
effect  on  the  final  results  for  a  wide  range  of  impact  velocities,  then 
Chat  parameter  should  be  eliminated  as  a  variable  from  the  model.  The 
two  proof  test  scenarios  described  in  Task  3  may  serve  as  a  good  vehicle 
for  these  parameter  variation  studies. 

Other  methods  involve  simplifying  Che  characterization  procedure. 
Using  the  shear  banding  damage  model  to  describe  the  behavior  of  a 
particular  material  requires  determining  which  parameters  govern  the 
kinetics  of  shear  banding  in  that  material.  Obtaining  these  parameters 
has  been  a  lengthy  procedure,  including,  in  addition  to  standard  quasi- 
static  tensile  tests  (if  unavailable  in  the  literature),  a  series  of  CKC 
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experiments  and  computer  simulations.  The  complete  procedure  has  been 
done  for  only  two  materials:  4340  steel  (Rc40)  and  an  RHA.  Computa¬ 
tional  simulation  of  scenarios  for  ocher  materials  would  be  much  easier 
if  the  characterization  procedure  were  simplified. 

Simplification  may  take  several  forms.  First,  the  number  of  CFC 
experiments  to  be  performed  may  be  reduced  by  modifying  the  experimental 
design.  In  the  past,  most  CFC  specimen  tubes  were  internally  loaded 
with  a  single  density  of  explosives,  and  the  density  was  varied  for 
different  experiments  to  obtain  a  wide  range  of  shear  band  damage. 
However,  in  one  experiment  performed  last  year  with  RHA,  the  explosive 
density  was  varied  as  a  function  of  axial  position  within  the  tube,  and 
an  adequately  wide  range  of  damage  was  obtained  from  that  single 
experiment. 

Second,  the  amount  of  data  to  be  obtained  from  a  series  of  CFC 
experiments  may  be  reduced.  If  the  parameter  variation  study,  outlined 
above,  has  shown  that  a  particular  material  parameter  has  a  negligible 
effect  on  the  simulated  outcome,  then  that  parameter  need  not  be 
measured.  For  example,  Che  shape  of  a  shear  band  in  its  slip  plane  is 
obtained  from  CFC  experiments,  and  in  fact,  a  significant  metallographic 
effort  is  needed  just  to  obtain  that  one  parameter.  In  the  three 
different  steels  studied  so  far,  the  shape  is  semicircular.  Substantial 
experimental  effort  can  be  saved  if  highly  eccentric  semi-elliptical 
shapes  would  produce  a  similar  outcome  in  the  computational  simulations. 

Also,  if  experimental  studies  of  a  class  of  materials  reveal  no 
significant  variation  in  a  particular  parameter,  then  that  parameter  may 
not  need  to  be  measured  for  a  new  material  of  that  class.  This  would  be 
true  even  if  varying  that  parameter  affected  the  computational  outcome. 
Using  the  same  example  as  above,  we  may  decide  that  all  steels  are 
likely  to  exhibit  semicircular  shear  bands.  There  is  of  course  a  risk, 
in  doing  this,  but  the  effort  saved  may  make  the  risk  worthwhile. 

Third,  we  may  not  need  to  follow  the  entire  procedure  for  charac¬ 
terizing  a  new  material  if  certain  of  its  properties  are  similar  to 
those  of  previously  characterized  materials.  This  simplification  could 


be  fundamentally  important.  If  we  could  identify  certain  material 
properties  from  which  we  could  (perhaps  in  conjunction  with  other  data) 
obtain  shear  band  kinetics  parameters,  then  we  would  not  only  have 
simplified  the  characterization  process,  but  also  increased  our  under¬ 
standing  of  the  shear  banding  phenomenon. 

For  example,  suppose  we  tested  a  series  of  CFC  specimen  tubes  that 
were  from  the  same  batch  of  steel  but  heat-treated  to  different  hard¬ 
nesses.  Suppose  also  that  the  results  indicated  that  certain  parameters 
such  as  shear  banding  nucleation  rate,  were  identical  over  the  range  of 
hardness  studied,  while  others,  such  as  nucleation  threshold  strain, 
varied  in  a  manner  that  could  be  related  to  other  material  variables 
(e.g.,  the  dynamic  yield  curve).  We  could  then  hypothesize  that,  at 
least  for  this  particular  steel,  the  nucleation  rate  was  an  invariant  of 
the  yield  curve,  while  the  threshold  strain  was  a  particular  function  of 
the  yield  curve.  And  the  hypothesis  could  be  tested  by  comparison  with 
CFC  results  of  yet  another  hardness  of  the  same  specimen  material. 
Finally,  if  similar  results  were  obtained  for  several  different  material 
within  a  particular  class  of  materials  (e.g.,  all  steels),  we  could  then 
generalize  the  results  to  all  materials  within  that  class  and  thus 
reduce  further  characterization  efforts. 

This  procedure  will  require  a  significant  effort  and  a  large  data 
base.  However,  we  do  have  an  appropriate  data  base  with  which  to  begin 
Che  procedure:  CFC  test  data  on  three  different  hardnesses  of  4340 
steel.  We  have  generated  parameters  for  one  of  the  hardnesses  (Rc40), 
but  have  yet  to  complete  data  analysis  for  hardnesses  Rc52  and  Rc2l. 

One  other  point  should  be  discussed  here:  using  alternative 
characterization  experiments.  Depending  upon  the  results  of  the 
simplification  efforts  described  above,  CFC  tests  may  not  be  the  most 
efficient  means  of  generating  shear  banding  parameters  for  certain 
materials.  For  example,  if  it  can  be  shown  that  for  a  particular  steel 
only  the  shear  band  threshhold  strain  and  the  dynamic  yield  curve  need 
to  be  determined,  using  an  alternative  experiment,  the  symmetric  rod 
impact  test,  would  reduce  the  total  effort.  This  latter  test  generates 


a  dynamic  yield  curve,  and  is  also  sensitive  to  the  shear  band  threshold 
strain. 


Task  5: 


jlication  of  Computational  Method  to  Various  Scenarios 


As  the  modif ications  described  in  the  first  four  tasks  are 
implemented  and  tested,  the  calculational  method  will  become  applicable 
to  an  ever  widening  variety  of  penetration  or  other  fragmentation 
scenarios.  Therefore,  we  recommend  appropriate  applications  of  the 
computational  method,  such  as  those  described  below. 


The  first  such  application  is  a  study  of  scaling  in  armor 
penetration.  Small-scale  testing  is  widely  used  by  armor  and  penetrator 
designers  because  it  is  time  and  cost  efficient.  It  is  thus  very  impor¬ 
tant  to  know  how  well  armor  penetration  and  the  resultant  fragmentation 
scale  with  size.  It  has  been  shown  previously  that  tensile  failure  does 

not  obey  replica  scaling  laws,  since  the  nucleation  and  growth  of  voids 

2 

are  rate-dependent  processes.  A  particular  scenario  that  results  in 
negligible  void  growth  for  a  small-scale  simulation  will  often  show 
substantial  void  growth  for  a  full-scale  simulation. 


Since  the  shear  band  nucleation  and  growth  criteria  in  SHEAR3  or 
SHEAR4  are  time-dependent  in  a  manner  similar  to  that  for  tensile  voids, 
it  is  expected  that  shear  banding  would  also  not  obey  replica  scaling 
laws.  A  preliminary  computational  study  of  scaling  in  CFC  experiments^1 
showed  a  significant  variation  of  shear  banding  for  different  scale 
factors. 


A  more  definitive  shear  band  scaling  study,  and  one  more  directly 
relavant  to  armor  penetration,  would  involve  selecting  a  set  of  pene¬ 
tration  scenarios  that  cover  a  wide  range  of  impact  velocities  and 
resultant  damage  (the  two  experiments  recommended  in  Task  3  for  the 
proof  tests  would  be  appropriate),  varying  the  scale  factor  in  both 
experimental  tests  (or  use  existing  data,  if  available)  and  computer 
simulations  of  these  scenarios,  and  observing  how  closely  the  simula¬ 
tions  predict  the  resultant  damage  for  different  scale  factors.  It  may 
turn  out,  for  example,  that  fragmentation  scales  reasonably  well  at  high 
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impact  velocities  but  that  damage  by  simple  plugging  at  lower  velocities 
does  not.  This  type  of  information  would  be  of  great  value  in  planning 
design  experiments. 


Other  possible  applications  include  the  following. 

DU-into-RHA  Penetration  Studies 

BRL  has  a  large  DU-into-RHA  fragmentation  data  base, 
covering  a  wide  range  of  velocity  and  obliquity, 
available  for  use  in  conjunction  with  computational 
studies.  The  simulation  of  previously  performed  DU  CFC 
experiments  must  be  completed  to  obtain  a  complete  set 
of  DU  shear  banding  parameters. 

Fragmentation  Device  Studies 

The  simulation  of  fragmentation  devices  has  been  limited 
by  the  fact  that  different  regions  fail  by  different 
mechanisms.  In  fragmenting  shells,  shear  banding 
predominates  near  the  interior  of  the  shell  wall,  while 
tensile  failure  predominates  at  the  exterior.  This 
limitation  will  be  greatly  diminished  once  the  feature 
(described  in  Task  l)  is  implemented  that  allows  a 
computational  cell  to  select  the  failure  appropriate  to 
its  loading  history. 


Task  6:  Documentation  and  Transfer  to  BRL 
<  Maximizing  the  potential  utility  of  a  computational  method  for 

i 

j  armor  penetration  can  be  achieved  only  by  transferring  to  BRL  the 

|  capability  to  use  the  method  routinely.  We  recommend  that  this  be  done 

by  thoroughly  documenting  the  C-HEMP  code  and  the  failure  models,  and  by 
directly  assisting  the  specific  BRL  personnel  who  will  be  responsible 
for  using  the  computational  method. 

Although  C-HEMP  was  developed  for  use  on  SRI's  VAX  computer,  only 
relatively  minor  modifications  will  be  needed  to  make  Che  code  compatible 
with  BRL's  CDC-6600  computer.  Depending  upon  the  final  size  of  C-HEMP, 
some  segmentation  may  be  required,  to  fit  the  memory  space  available  on 
the  BRL  computer.  Also,  some  FOPTRAN  statements  will  need  to  be  changed 
for  language  compatibility. 
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